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1 Introduction 

One of the oldest but still actual problems in nuclear physics is related to the determination of the 
interaction between the nucleons. A quantitative understanding of the nuclear force is crucial in order 
to describe the properties of nuclei and nuclear matter in terms of hadronic degrees of freedom. The 
conventional way to parametrize the nuclear force utilizes the meson-exchange picture, which goes back 
to the seminal work by Yukawa [1]. His idea, followed by the experimental discovery of tt- and heavier 
mesons (p, u, . . . ), stimulated the development of boson-exchange models which still provide a basis 
for many modern, highly sophisticated phenomenological nucleon-nucleon (NN) potentials. 

According to our present understanding, the nuclear force is due to residual strong interactions between 
the color-charge neutral hadrons. A direct derivation of the nuclear force from QCD, the underlying 
theory of strong interactions, is not yet possible, see however Ref. [2] for a recent attempt using lattice 
QCD. In order to provide reliable input for few- and many-body calculations, a (semi-) phenomenological 
approach has been followed over the past few decades aiming to achieve the best possible description of 
the available low-energy NN data. As will be discussed in section [2j the two-nucleon potential can be 
decomposed in only few different spin-space structures, so that the corresponding radial functions can 
be parameterized using an extensive set of data. Although the resulting models provide an excellent 
description of experimental data in many cases, there are certain major conceptual deficiencies that 
cannot be overcome. In particular, one important concern is related to the problem of the construc- 
tion of consistent many-body forces. These can only be meaningfully defined in a consistent scheme 
with a given two-nucleon interaction j3]. Notice that because of the large variety of different possible 
structures in the three- nucleon force, following the same phenomenological path as in the NN system 
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and parametrizing its most general structure seems not to be feasible without additional theoretical 
guidance. Clearly, the same problem of consistency arises in the context of reactions with electroweak 
probes, whose description requires the knowledge of the corresponding consistent nuclear current op- 
erator. Further, one lacks within phenomenological treatments a method of systematically improving 
the theory of the nuclear force in terms of the dominant dynamical contributions. Finally, and most 
important, the phenomenological scheme provides only a loose connection to QCD. 

Chiral perturbation theory (ChPT) is an effective field theory (EFT) of QCD which exploits its sym- 
metries and symmetry-breaking pattern and allows to analyze the properties of hadronic systems at low 
energies in a systematic and model independent way. We will see in section [3] that QCD with two fla- 
vors of the u- and d-quarks and, to a less extent, with three flavors of the u-, d- and s-quarks, exhibits 
an approximate chiral symmetry which is explicitly broken due to non-vanishing (but small) quark 
masses. In addition, the chiral symmetry is also spontaneously broken down to its vector subgroup. 
These symmetry/symmetry- breaking pattern manifest themselves in the hadron spectrum leading, in 
particular, to a natural explanation of the very small masses (compared to other hadrons) of pions which 
play the role of the corresponding Goldstone bosons. Most important, the Goldstone boson nature of 
the pions implies that they interact weakly at low energy and allows to calculate low-energy observables 
in the pion and single-nucleon sector in perturbation theory. The situation in the few-nucleon sector 
is conceptually much more complicated due to the strong nature of the nuclear force which manifests 
itself in the appearance of self-bound atomic nuclei and invalidates a naive application of perturbation 
theory. As pointed out by Weinberg, the breakdown of perturbation theory in the few-nucleon sector 
can be traced back to the infrared enhancement of reducible time-ordered diagrams which involve purely 
nucleonic intermediate states and can be resummed by iterating the corresponding dynamical equation 
[HE]. These important observation made in Weinberg's seminal papers opened a new era in nuclear 
physics and has triggered an intense research activity along these lines. In these lectures I will outline 
the basic concepts of chiral effective field theory and its application to nucleon-nucleon scattering and 
the derivation of the nuclear force. 

The manuscript is organized as follows. In sectional discuss the general structure of the nuclear force 
and outline the main ingredients of the conventional NN potentials. Section [3] provides an elementary 
introduction to chiral perturbation theory. Generalization of EFT to strongly interacting nuclear sys- 
tems is discussed in section |4j Derivation of the nuclear forces in chiral EFT is outlined in section [H 
A brief summary is given in section [6j 

2 Nuclear potentials and nucleon-nucleon scattering 

The most general structure of a non-relativistic two-nucleon potential is expressible in terms of just 
a few operators. The potential can be viewed as an operator acting in the position, spin and isospin 
spaces of the nucleons. It is instructive to discuss its isospin structure separately from the operators 
acting in the position-spin space. 

The isospin structure of the two-nucleon force falls into the four different classes according to the 
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classification of Ref. [6] : 



Class I: 



V\ = ai + A ti ■ t 2 , 
Vn = a n if r| , 

Vm = "in (ti + t|) , 

Viv = aiv (tx 3 - rf ) + fty [n x r 2 ] 3 . 



Class II: 



Class III: 



(2.1) 



Class IV: 



Here, «i, /3i are position-spin operators and Ti are Pauli isospin matrices of a nucleon i. The operator 
/3iv has to be odd under a time reversal transformation. While class (I) forces are isospin-invariant, all 
other classes (II), (III) and (IV) are isospin-breaking. Class (II) forces, Vn, maintain charge symmetry 
but break charge independence. They are usually referred to as charge independence breaking (CIB) 
forces. Charge symmetry represents invariance under reflection about the 1-2 plane in charge space. 
The charge symmetry operator P cs transforms proton and neutron states into each other and is given 
by P cs = e iwT2 with T = J2i T i/ 2 bein g the total 

isospin operator. Class (III) forces break charge 
symmetry but do not lead to isospin mixing in the NN system, i.e. they do not give rise to transitions 
between isospin-singlet and isospin-triplet two- nucleon states. Finally, class (IV) forces break charge 
symmetry and cause isospin mixing in the NN system. 

Exercise: show that class-Ill two-nucleon forces do not lead to isospin mixing in the two-nucleon 
system, i.e. they commute with the operator T 2 . Does this still hold true for systems with three and 
more nucleons? 

Let us now discuss the position-spin structure of the potential. For the sake of simplicity, I restrict 
myself to the isospin-invariant case. The available vectors are given by the position, momentum and spin 
operators for individual nucleons: fi, f 2 , pi, p 2 , <?i, <?2- The translational and Galilean invariance of the 
potential implies that it may only depend on the relative distance between the nucleons, r = f\ — r 2 , 
and the relative momentum, p = (p\ — p 2 )/2. Further constraints due to (i) rotational invariance, 
(ii) invariance under a parity operation, (hi) time reversal invariance, (iv) hermiticity as well as (v) 
invariance with respect to interchanging the nucleon labels, 1 2, lead to the following operator form 
of the potential [7] : 



Throughout this work, two-nucleon observables will be computed by solving the Lippmann-Schwinger 
equation in momentum space. It is, therefore, instructive to look at the momentum-space representation 
of the potential, V(p', p) = (p'\ V\p), with p and p' denoting the two-nucleon center of mass momenta 
before and after the interaction takes place. Following the same logic as above, the most general form 
of the potential potential in momentum space can be shown to be: 




where L = f x p, S = {5\ + <?2)/2 and S\2(x ) = 3<ti ■ x u% • x — <j\ • ai with x = x/\x\. The operators 
entering the above equation are multiplied by scalar operator- like functions that depend on r 2 , p 2 and 



L 2 . 



spim 



<?i • &2, Si 2 (q), Si 2 {k), iS ■ qxk, o\ ■ qx ka 2 ■ qx k 



}x{l 




(2.3) 
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where q = p 1 
depend on p 2 

that enter Eq. (2.3) denote the corresponding eigenvalues. It should also be emphasized that further 



— p and k = (p 1 +p)/2. The operators a re m ultiplied with the scalar functions that 
p' 2 and p ■ p' . Notice that contrary to Eq. (2.2) which involves the operator p, p and p' 



spin- momentum operators contribute in the case of class- IV isospin-breaking interactions. 

For low-energy processes I will be focused in here, it is convenient to switch to the partial wave basis 
\p) — > \plmi). A two-nucleon state \p(ls)jm,j) in the partial- wave basis depends on the orbital angular 
momentum I, spin s, the total angular momentum j and the corresponding magnetic quantum number 



rrij. The partial wave decomposition of the potential in Eq. (2.3) is given by 



with 



(p' (I' s')j'm' j \V\p(ls)jm j ) 



Ss'sVyHp', p) , 



(2.4) 



V$(p',p) = ^2 I d P' dpc(l',s,j;mi,mj - m'^rrij) c(l,s,j;mi,mj -m h mj) 

m! , mi 

x Y Vm'(P') YimXP) ( sm j ~ m'^Vip', p)\smj - mi) , 



(2.5) 



where c(l,s,j;mi,mj — mi,mj) are Clebsch-Gordan coefficients and Yi mi (p) denote the spherical har- 



monics. The first two Kronecker <5's on the right-hand side of the first line in Eq. (2.4) reflect the 



conservation of the total angular momentum. Rotational invariance of the potential prevents the de- 
pendence of the matrix elements on the magnetic quantum number rrij . The conservation of the total 



spin of the nucleons can be easily verified explicitly for all operators entering Eq. (2.3). I stress, how 



ever, that transitions between the spin-singlet and spin-triplet channels are possible in a more general 



case of the broken isospin symmetry. For each individual operator entering Eq. (2.3), the expression 



(2.5) can be simplified and finally expressed as an integral over p-p' with the integrand being written in 



terms of the corresponding scalar function and Legendre polynomials. Explicit formulae can be found 
e.g. in [8], see also Ref. [9] for a recent work on this topic. 

The Lippmann-Schwinger (LS) equation for the half-shell T-matrix in the partial wave basis has the 
form 



T»(p',p) = v;i(p', P ) + 



E 

I" 



00 dp"p" 2 



(2tt) e 



Vnl,{p',p") 



m 



p + if] 



(2.6) 



with m denoting the nucleon mass and r\ — > + . In the uncoupled case, I is conserved. The relation 
between the on-shell S- and T-matrices is given by 



s vM = § i'i ~ ^2P mT n (?) • 
The phase shifts in the uncoupled cases can be obtained from the 5-matrix via 



(2.7) 



cxp 



2^ 



S )] =exp(2^ 



f2i 
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where I use the notation Sf J . The so-called Stapp parametrization of the S'-matrix in the coupled 
channels (j > 0) is defined as: 



S 




cos (2e) exp (2«<5 J -f_ 1 ) i sin (2e) exp (iSjLi + ^j+i/ 
i sin (2e) exp (i^i + ^j+i) cos (2e) exp (2i<5^ 1 ) 



and is related to another frequently used parametrization due to Blatt and Biedenharn in terms of 5 
and e via the following equations: 

? ? • r? r \ tan(2e) . ,~ ? . sin(2e) n . 

^•-1 + ^+1 = ^-1 + ^+1, sm(^_i - <y i+ i) = tan 2g ( , sm^.i - = — — . (2.9) 

The appearance of the electromagnetic interaction requires special care when calculating scattering 
observables due to its long-range nature. In particular, the S-matrix has to be formulated in terms 
of asymptotic Coulomb states. The electromagnetic interaction between the nucleons is driven by the 
Coulomb force and, to a lesser extent, magnetic moment interactions and vacuum polarization. It 
should also be emphasized that the expansion of the scattering amplitude in partial waves converges 
very slowly in the presence of the magnetic moment interactions. For explicit expressions and a detailed 
discussion on their implementation when calculating nucleon-nucleon observables the reader is referred 

to nm. 



The deuteron wave function and binding energy Ed are obtained from the homogeneous part of Eq. ( 2.6 ): 

where s = j = 1, I = I' = 0, 2. Once phase shifts are calculated, nucleon-nucleon scattering observables 
can be computed in a standard way, see [ITJ [12] . 

The appearance of only a few structures in the most general expression for the two-nucleon force, see 



Eq. (2.3), and the large amount of available low-energy nucleon-nucleon scattering data motivated and 
enabled the development of modern high-precision phenomenological potential models such as e.g. the 
CD-Bonn 2000 pj, Argonne Vis (AV18) [Hj and Nijmegen I, II potentials [EJ. The general strat- 
egy involves incorporating the proper long-range behavior due to the electromagnetic interaction and 
the one-pion exchange potential which is important to correctly describe the low-energy behavior of 
the amplitude, cf. section [4~2| and parametrizing the medium- and short-range contributions in a gen- 
eral way. AV18, a local r-space potential, can be viewed as a representative example. It includes (i) 
electromagnetic interactions multiplied by short-range functions to account for the finite size of the 
nucleon, (ii) regularized one-pion exchange potential including isospin-breaking corrections due to dif- 
ferent masses of the charged and neutral pions, (iii) some additional phenomenological isospin-breaking 
terms of a shorter range, (iv) medium-range (short-range) contributions of Yukawa-type (Woods-Saxon 
type) multiplying the operators in Eq. ( |2.2[ ). With about 40 adjustable parameters, it describes the 
proton-proton and neutron-proton scattering data with Xdatum = 1-09- Other high-precision potentials 
are constructed in a similar way and allow to reproduce the data or phase shifts from e.g. the Nijmegen 
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Tiab(MeV) T, a0 (MeV) T, aa (MeV) 



Figure 1: 3 Si (left panel) and 3 D\ (right panel) phase shifts and the mixing angle ei (middle panel) 
calculated from several modern high-precision potentials in comparison with the results of the Nijmegen 
PWA. The phase shifts and the mixing angle are shown in degrees. Plots are generated through the 
NN-Online web site http://nn-online.org 



partial wave analysis (PWA) with a comparable accuracy. This is visualized in FigjT] I refer the reader 
to Ref. |16j for a recent review article on the modern high-precision potentials. 

While various phenomenological potentials provide an accurate representation of the nucleon-nucleon 
phase shifts and most of the deuteron properties, the situation is much less satisfactory when it comes to 
the much weaker but necessary three-nucleon forces. Such three-body forces are needed to describe the 
nuclear binding energies and levels, as most systematically shown by the Urbana-Argonne group [17J. 
Systematic studies of the dynamics and reactions of systems with three or four-nucleons further sharpen 
the case for the necessity of including three-nucleon forces, see e.g. |18j . A phenomenological path to 
modeling the three-nucleon force following the same strategy as in the two-nucleon case seems to be 
not feasible (at least, at present). Indeed, in the case of two nucleons, the potential can be decomposed 
in only a few different spin-space structures, and the corresponding radial functions can be adjusted to 
the extensive set of data. Such an approach would, however, fail for the three-nucleon force due to the 
large variety of different possible structures, a scarcer data base and considerably more time consuming 
calculations required. 

While the conventional approach based on the high-precision two-nucleon potentials accompanied with 
the existing three-nucleon force models enjoyed many successes and is frequently used in e.g. nuclear 
structure and reaction calculations, it remains incomplete as there are certain deficiencies that can only 
be overcome based on EFT approaches. These are: (i) it is very difficult - if not impossible - to assign a 
trustworthy theoretical error, (ii) gauge and chiral symmetries are difficult to implement, (iii) none of the 
three-nucleon forces is consistent with the underlying nucleon-nucleon interaction models/approaches 
and (iv) the connection to QCD is not at all obvious. 
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Figure 2: A localized charge distribution generates an electrostatic potential which can be described 
in terms of the multipole expansion. 

3 Chiral perturbation theory: An elementary introduction 

Effective field theories have proved to be an important and very useful tool in nuclear and particle 
physics. One understands under an effective (field) theory an approximate theory whose scope is to 
describe phenomena which occur at a chosen length (or energy) range. The main idea of this method 
can be illustrated with the following example from classical electrodynamics. Consider a localized 
charge distribution in space of a size a. The resulting electrostatic potential at any given position R 
can be calculated by integrating over the elementary charges and using the familiar expression for the 
Coulomb potential generated by a point charge: 







(3.12) 



with i, j denoting the Cartesian 



components allows to rewrite the integral as 




pif) 



I? + ~r3 X] RiPi + am ^(ZRiRj ~ SijR 2 )Qij + 



(3.13) 



R-f 



i ij 
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where the total charge q, dipole moment Pi and quadrupole moment Qij are defined via 

d 3 rp(r), Pi= f d 3 rp(f)n, Q fj = [ d 3 r p(f )(3r i r j - 5 ijr 2 ) . (3.14) 



The expression in Eq. (3.13) represents the well-known multipole expansion for the electrostatic po- 



tential. When truncated, it provides an approximation to the "underlying theory" given by the exact 



expression (3.11). The multipoles entering every term in this expansion contain certain amount of 
information about the charge distribution and can, of course, be calculated provided p{r) is known. 
The multipole expansion is, however, particularly useful if p{r) is unknown (except for the fact that 
it is localized). It then allows to describe the electrostatic potential at every point in space far from 
the charge distribution with, in principle, an arbitrarily high accuracy provided one has enough data 
(e.g. experimentally measured values of the electrostatic potential at some points) to determine the 
desired number of the multipoles. 

Chiral Perturbation Theory (CHPT) is the effective theory of QCD (more generally, of the Standard 
Model) which was formulated by Weinberg |19] and developed in to a systematic tool for analyzing low- 
energy hadronic observables by Gasser and Leutwyler |20^ 121] . In this section, I give a brief overview 
of the foundations of this approach. My main purpose here is to outline the logical steps which are 
needed in order to set up this theoretical framework. I will also give references to the existing extensive 
literature on this subject which is suitable for further reading. 



3.1 Chiral symmetry of QCD 

Symmetries provide powerful constraints on effective interactions and thus play the crucial role for 
effective field theories. In the following, I will discuss the symmetries of QCD which are relevant in the 
context of ChPT. Consider the QCD Lagrangian in the two-flavor case of the light up and down quarks 

£ Q CD = q - M)q - ~G%G a ^ , (3.15) 

where = — ig s G®T a with T a , (with a = 1 ... 8) are the SU(3) co i or Gell-Mann matrices and q the 
quark fields. Further, G^ u are the gluon field strength tensors, and the quark mass matrix is given by 
M. = diag(m M , m^). I do not show in Eq. (3.15) the 6- and gauge fixing terms which are not relevant for 
our consideration. It is instructive to write the QCD Lagrangian in terms of the left- and right-handed 
quark field components defined by q R = (1/2)(1 + 75)9, qi = (1/2)(1 — 75)9: 

£ Q CD = q L ipQL + qRilpqR ~ q L Mq R - q R Mq L - ^G^G^ . (3.16) 

We see that the left- and right-handed quark fields are only connected through the mass term. Given 
the smallness of the light quark masses |22j Q 

m u ~ 1.5 . . . 3.3 MeV, m d ~ 3.5 . . . 6.0 MeV, (3.17) 



1 The following values correspond to the MS scheme at scale fj, = 2 GeV. 
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as compared to the typical hadron masses of the order of 1 GeV, the quark mass term can, to a good 



approximation, be neglected. The Lagrangian in Eq. (3.16) is, therefore, approximately invariant under 



independent global flavor rotations of the left- and right-handed quark fields: 



<1L 



Lq L = exp {-iO L ■ r/2)q L , 



q' R = Rq R = exp {-iO R ■ r/2)q R , 



(3.18) 



where r denote the Pauli matrices in the flavor space and Ql,r are the corresponding rotation an- 
gles. The corresponding symmetry group SU(2)l x SU(2) r is referred to as the SU(2) chiral group. 
According to Noether's theorem, there are six conserved currents 



L] 



QRln^QR 



(3.19) 



which can equally well be expressed in terms of the vector and axial- vector currents V* 
A 1 = Rt — L % . The corresponding conserved charges generate the algebra of the chiral group 



L). + R). and 



Q}, Q\ = ie ijk Q k j with I = L,R, 



Ql-i Or 



or, equivalently, 



Qyi Qy 



ie ijk Q k v , Q% Q\ = ie ijk Q k v , Q v , Q\ = ie ijk Q k - 



l A ■ 



(3.20) 



(3.21) 



Application of the above commutation relations to hadronic reactions was at the heart of the current 
algebra calculations in the early seventies of the last century. 

The Lagrangian for massless u- and d-quarks is, in fact, invariant under even a larger group of trans- 
formations in the flavor space, namely SU(2)i x SU(2) R x U(l)y x U(1)a- While the vector U(l) 
corresponds to quark number conservation, the axial U(1)a current is known to be broken by quantum 
effects (the so-called U(1)a anomaly) and thus does not represent a symmetry of the quantum theory. 

In spite of the fact that QCD for two light flavors is approximately chiral invariant, its ground state 
is not symmetric with respect to SU(2)l x SU(2) r but only with respect to its vector subgroup 
SU(2)v C SU(2)i x SU(2)r generated by the charges {Qy}. This means that the axial charges 
do not annihilate the vacuum, that is Qy|0) = while Q^|0) ^ 0. Evidence of the spontaneous 
breakdown of the chiral symmetry comes from various sources. For example, hadrons occur in nearly 
degenerate isospin multiplets corresponding to SU(2)y which implies that this group is realized in the 
usual Wigner-Weyl mode. If this were the case for the chiral group, one would observe larger chiral 
multiplets containing particles of opposite parity since the charges Q l v and Q A have opposite parity. 
Generally, no such parity doubling is observed in the hadron spectrum. Another strong argument in 
favor of the spontaneous breakdown of the chiral symmetry comes from the existence of unnaturally 
light (in comparison with other hadrons) pseudoscalar mesons (pions) being natural candidates for the 
corresponding Nambu-Goldstone bosons. Pions are not exactly massless but acquire a small mass due to 
the explicit chiral symmetry breaking by the nonvanishing quark masses. These and further arguments 
coming from both the theory and experiment indicate undoubtly that the chiral SU(2)l x SU(2) r 
group is spontaneously broken down to SU(2)y. 
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I now pause to summarize the content of this section. QCD Lagrangian in the two-flavor case of the 
up- and down-quarks is approximately invariant under global chiral SU(2)l x SU(2)r transformations. 
The chiral symmetry is broken explicitly due to the nonvanishing quark masses, m u 7^ 0, rrid 7^ 0. In 
addition to this explicit symmetry breaking, SU (2)l x SU(2)r is also broken spontaneously down to 
the isospin group SU(2)y. The three corresponding pseudoscalar Goldstone bosons are identified with 
pions whose small masses emerge due to nonvanishing quark masses. There exist a large mass gap in 
the hadron spectrum: M p ~ 770 MeV »M T ~ 140 MeV. 



3.2 Effective Lagrangian for Goldstone bosons 

We now turn to the effective description of low-energy QCD dynamics. The simplest possible case 
emerges when the energy is chosen so small that only pions need to be treated as explicit degrees 
of freedom. All other hadrons are much heavier and can be integrated out from the theory. The 
main ingredient of any effective field theory is the most general effective Lagrangian that involves all 
possible terms which are consistent with the symmetries of the underlying theory. Let us, for the time 
being, consider the so-called chiral limit of QCD, i.e. the idealized world in which quarks are massless 
and the chiral symmetry of Cqcd is exact. The task is then to construct the most general chiral 
invariant Lagrangian for pion fields. In order to do that we first need to figure out how pions transform 
with respect to chiral rotations. Our knowledge of the pion transformation properties with respect to 
SU(2)l x SU(2)r can be summarized by the following two observations: 

• Pions build an isospin triplet and thus transform linearly under SU(2)y C SU(2)i x SU(2)r 
according to the corresponding irreducible representation; 

• The chiral group must be realized nonlinearly. This follows immediately from the geometrical 



argument based on the fact that the Lie algebra of SU (2) l x SU (2)r in Eq. (3.20 ) is isomorphic to 
that of SO (4). We know that one needs three coordinates in order to construct the smallest non- 
trivial representation, the so-called fundamental representation, of the three-dimensional rotation 
group. Similarly, the smallest nontrivial representation of the four-dimensional rotation group 
S"0(4) is four-dimensional. We have, however, only three "coordinates" at our disposal (the 
triplet of the pion fields) ! 

To construct a non- linear realization of SO (A) we begin with the usual representation describing four- 
dimensional rotations of a vector (tt, a) = (ttj, vr2, tt^, a). For an infinitesimal rotation parametrized 
by six angles {9^' A }, with i = 1,2, 3, we have: 



7T \ 50(4) 

a ) \ <r 



71 1 (3.22) 
a ' 
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where 



E^ 



i=l 
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(3.23) 



0$ / 



Notice that the set of rotations generated by V% builds a subgroup of 50(4), namely the group of three- 
dimensional rotations 50(3) C 50(4) which is locally isomorphic to SU{2). The four real quantities 
(tv, a) define the smallest nontrivial chiral multiplet and represent the field content of the well-known 
linear sigma model. To switch from the above linear realization (i.e. representation) of 50(4) to the 
nonlinear one, we observe that, in fact, only three of the four components of (tv, a) are independent with 
respect to four-dimensional rotations. These three independent components correspond to coordinates 
on a four-dimensional sphere since tv and a are subject to the constraint 



2 i 2 
TV + a 



(3.24) 



where F is a constant of dimension mass. Making use of this equation to eliminate a in Eq. (3.22) we 
end up with the following transformation properties of tv under 50(4): 



TV 



TV 



TV 1 = TV + X TV . 



TV' = TV + 9 A ^F 2 -TV 2 , 



(3.25) 



where 6 V ' A = {@Y ,A } with i = 1, 2, 3. The nonlinear terms (in tv) on the right-hand side of the second 
equation give rise to the nonlinear realization of 50(4). This is exactly what we wanted to achieve: the 
chiral group SU(2)l x SU(2)r ~ 50(4) is realized nonlinearly on the triplet of pions which, however, 
transform linearly under isospin SU(2)y ~ 50(3) rotations parametrized through the angles {Oy}. 

As a last remark note that the four-dimensional rotations of (tv, a) can be conveniently written using 
the 2x2 matrix notation by introducing the unitary matri^ 



U = — (crl 2x2 +itv -t) , 
and demanding the transformation properties of U under chiral rotations to be: 

U — > U' = LUR ] . 
Here, L and R are SU(2)l x SU(2)r matrices defined in Eq. ( |3.18 ). 



(3.26) 



(3.27) 



Exercise: verify that infinitesimal transformations of (tv, a) induced by Eq. (3.21) with 
O v = (9ft + Ol)/2 and A = (Or — Ol)/2 have indeed the same form as the ones given in 



Eq. (3.22). 



For U to be unitary, a and tv have to fulfill Eq. (3.241 
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Clearly, the transition to the nonlinear realization is achieved by 



U = - (a l 2x2 + iw ■ t) — ► U=j (Vi ?2 -7r 2 l2x2 + t7T • t) , (3.28) 

leaving pions as the only remaining degrees of freedom. Notice that the ground state of the theory 
is characterized by a vanishing vacuum expectation values of 7r and corresponds to a particular point 
on the considered four-dimensional sphere (one of two crossing points between the sphere and the 
ex-axis). In accordance with the spontaneous breaking of chiral symmetry, it not SO (4)- but only 
50(3) ~ 5C/(2) y -invariant. 

It is now a simple exercise to construct the most general chiral-invariant Lagrangian for pions in terms 
of the matrix U. The building blocks are given by U, and derivatives of these quantities. Notice 
that since I consider here only global chiral rotations, i.e. L and R do not depend on space-time, the 



quantities like e.g. d^d u U transform in the same way as U itself, i.e. according to Eq. (3.27). Chiral 
invariant terms in the effective Lagrangian terms can be constructed by taking a trace over products of 
U, W and their derivatives. Lorentz invariance implies that the number of derivatives must be even, 
so that the effective Lagrangian can be written as 

£ 7r = 4 2 )+4 4 ) + ... . (3.29) 

Notice that is simply a constant since UU^ = l 2 x2- The lowest-order Lagrangian involves just a 
single term 

4 2 ' = j(a^[/t), (3.30) 

where (. . .) denotes the trace in the flavor space. Terms involving d^d^U or d^d^U^ are not independent 



and can be brought to the form of Eq. (3.30) by using partial integration. The constant F /4 ensures 



that the Lagrangian has a proper dimension (both F and ir have a dimension of mass) and is chosen 



in such a way that Eq. (3.30) matches the usual free Lagrangian for massless scalar field when written 
in terms of pions: 

4 2 ) = l -d,-K ■ 3"ir + ^ (d,n ■ tt) 2 + 0(tt 6 ) . (3.31) 
The a-priori unknown constants that accompany the terms in the effective Lagrangian are commonly 



called the low-energy constants (LECs), cf. the multipoles in Eq. (3.13). The LEC F can be identified 
with the pion decay constant in our idealized world (with quark masses being set to zero) [20J. In 
the real world, it is measured to be F^ = 92.4 MeV. Notice further that the pions are massless in the 
idealized world. Higher-order Lagrangians £-ir\ • • •, can be constructed along the same lines by 
using partial integration and equations of motion to eliminate the redundant terms. 

Let us now pause to summarize what has been achieved so far. We have explicitly constructed a 
particular nonlinear realization of SU(2)l x SU{2)r in terms of pion fields which, as desired, build a 
representation of SU (2)y C SU (2) i x SU (2)r and learned how to write down the most general possible 
effective Lagrangian. The constructed nonlinear realization of the chiral group is, however, not unique. 



Different realizations emerge by choosing different parametrizations of the matrix U in Eq. (3.27) in 
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terms of pion fields such as, for example, the exponential parametrization U = exp(i7r -t/F). Generally, 
only the first three terms in the expansion of U in powers of the pion field are fixed by unitarity, 

2 2 4 

U(tt) = l 2x2 + - ^ - + (8« - 1)^ + 0(tt 5 ) . (3.32) 

Here, a is an arbitrary constant which reflects the freedom in parametrizing the matrix U. This raises 
the concern that observables that one intends to compute in the EFT are possibly affected by this 
non-uniqueness which would be a disaster. Fortunately, this is not the case. As shown by Coleman, 
Callan, Wess and Zumino |23l I24j . all realizations of the chiral group are equivalent to each other 
modulo nonlinear field redefinitions 

vfi TT- = TTi F [it] with F[0] = 1. (3.33) 

According to Haag's theorem [25j, such nonlinear field redefinitions do not affect S-matrix elements. 

So far, I considered the chiral limit corresponding to the idealized world with the masses of the up- 
and down-quarks being set to zero. This is fine as a first approximation but, of course, one would like 
to systematically improve on it by taking into account corrections due to nonvanishing quark masses. 
For that, we have to include in the effective Lagrangian all possible terms which break chiral symmetry 
in exactly the same way as does the quark mass term in Cqcd- Consider, for example, the quark 
mass term with m u = rrid = m q ^ which breaks chiral but preserves isospin symmetry. Recalling 
the geometrical interpretation with the four-dimensional rotation group and coordinates (7r, a), the 
quark mass term can be viewed as a vector that points along the (0, cr)-direction. Its effects can be 
systematically taken into account by including in the effective Lagrangian not only S'0(4)-scalars but 
also the corresponding components of all possible S'0(4)-tensors and multiplying the resulting terms 
by the appropriate powers of m q , see [26] for more details and the explicit construction along these 
lines. A simpler (but equivalent) method makes use of the following trick. Consider the massless QCD 
Lagrangian in the presence of an external hermitian scalar field s interacting with the quarks via the 
term —qsq. The resulting Lagrangian is chiral invariant provided the scalar source s transforms under 
chiral rotations according to: 

s -> s' = LsR j = Rstf , (3.34) 

where the second equality follows from the hermiticity of s. To recover QCD from the new theory, 
the external field needs to be set to the value s = M. To account for the explicit chiral symmetry 
breaking, we first write down the effective Lagrangian for the new theory by listing all possible chiral 
invariant terms constructed from s and (derivatives of) U and W and then set s = M. Since the 
quark masses are treated as a small perturbation, the leading symmetry-breaking terms should contain 
a minimal possible number of derivatives and just one insertion of s. Given that U — > W under parity 
transformation, there exist only one symmetry-breaking term without derivatives: 



B 

= F 2 B(m u + m d ) - -{m u + m d )-K 2 + 0(n 4 ) , (3.35) 

s=M 1 



where B is a LEG The first term is a constant and does not contribute to the .S-matrix. The second 
one gives rise to the pion mass — (l/2)M 2 7r 2 with M 2 = {m u + rrid)B. Note that to leading order 
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in m Ut d, one has equal masses for all pions tt + , tt and ir . Further, the LEC B can be shown to be 
related to the quark condensate according to (0|uu|0) = (0|dd|0) = -F 2 B(1 + 0(M)) [20]. Modulo 
corrections of higher order in the quark masses, the experimentally measured pion mass M n coincides 
with M: M 2 = M 2 + 0{m 2 q ). 

How important are chiral-symmetry breaking terms as compared to chiral-invariant ones? When con- 
structing the most general chiral-invariant effective Lagrangian, various terms were classified according 
to the number of derivatives. When calculating observables, the derivatives acting on pion fields gener- 
ate powers of external momenta which are assumed to be low. The EFT considered so far is typically 
applied to processes characterized by pion momenta of the order of the pion massj^] It is, therefore, 
natural to count the pion mass in the effective Lagrangian on the same footing as a derivative. We thus 
end up with the following lowest-order Lagrangian: 

F 2 

4 2 ) = (dpUWrt + 2B{MU + MU*)) • (3.36) 
For the sake of completeness, the next-higher order Lagrangian reads [20 : 

4 4) = l j{d^ud^) 2 + l j(d^ud u u^(d^ud u u^) + ^(2BM(u + u^) 2 + ... 

- h-(2BM{U-U^)) 2 , (3.37) 
lb 

where the ellipses refer to terms that involve external sources and the Zj are the corresponding LECs. 



3.3 Power counting 

Having constructed the most general effective Lagrangian for pions in harmony with the chiral symmetry 
of QCD, we now need to figure out how to compute observables. At first sight, the effective Lagrangian 
seems to be of less practical value due to an infinite number of the unknown LECs. Even worse, all 
interaction terms entering C n are non-renormalizable in the usual sens^] contrary to field theories such 
as e.g. QED and QCD. What at first sight appears to cause a problem, namely non-renormalizability of 
the theory, in fact, turns out to be a crucial feature for the whole approach to be useful. As demonstrated 
in the seminal paper by Weinberg [19] , the effective Lagrangian C T can be used to compute low-energy 
observables (such as e.g. scattering amplitudes) in a systematically improvable way via an expansion 
in powers of Q/A x , where Q represents the soft scale associated with external momenta or the pion 
mass M n while Ay, the so-called chiral-symmetry-breaking scale, is the hard scale that drives the 
LECs in C-k. This expansion is referred to as the chiral expansion, and the whole approach carries 
the name of chiral perturbation theory. Consider an arbitrary multi-pion scattering process with all 
initial and final pion momenta of the order of M n . In order to decide on the importance of a particular 
Feynman diagram, we have to determine the power of the soft scale associated with it. For that we 

3 For an example of EFT in a different kinematical regime with nonrelativistic mesons see [27J . 

4 This implies that the structure of local ultraviolet divergences generated by loop diagrams with vertices from is 
different from . 
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Figure 3: Tadpole contribution from (left) and tree contribution from d£ (right) to the pion self 
energy. 



first need to clarify an important issue related to the counting of virtual momenta which are being 
integrated over in the loop integrals. What scale do we associate with such virtual momenta? When 
calculating Feynman diagrams in ChPT, one generally encounters two kinds of loop integrals. First, 
there are cases in which the integrand dies out fast enough when the loop momenta go to infinity so 
that the corresponding integrals are well-defined. Since the hard scale only enters the LECs and thus 
factorizes out, the integrands involve only soft scales (external momenta and the pion mass) and the 
loop momenta. Given that the integration is carried over the whole range of momenta, the resulting 
mass dimension of the integral is obviously driven by the soft scales. Thus, in this case we can safely 
count all virtual momenta as the soft scale. The second kind of integrals involves ultraviolet divergences 
and requires regularization and renormalization. Choosing renormalization conditions in a suitable way, 
one can ensure that virtual momenta are (effectively) of the order of the soft scale. This is achieved 
automatically if one uses a mass-independent regularization such as e.g. dimensional regularization 
(DR). Consider, for example, the integral 

/ = / (27T)4 /2 _ ^2 + i£ ' ( 3 - 38 ) 

that enters the pion self energy due to the tadpole diagram shown in Fig.[3j Evaluating this quadratically 
divergent integral in dimensional regularization one obtains 

1 - '" s - / WfW^m = t£ (f ) + 2M2i <"> + °< d " 4 » • < 3 - 39 > 

where je is the Euler constant and \x is the scale introduced by dimensional regularization and the 
quantity L{p) is given by 

m = iS { d^i - \ {H ^ ] + + 1} } ' = -°- 5772is ■ • • • ( 3 - 4 °) 

The second term on the right-hand side of the above expression diverges in the limit d — >■ 4 but can 
be absorbed into an appropriate redefinition of the LECs in (renormalization). Notice that, as 
desired, the mass dimension of the finite term is driven by the soft scale M. I further emphasize 
that the scale \i introduced by dimensional regularization has to be chosen of the order [i ~ M v in 
order to prevent the appearance of large logarithms in DR expressions. Last but not least, note that 
one can, in principle, use different regularization methods such as e.g. cutoff regularization provided 
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Figure 4: Diagrams contributing to pion-pion scattering in at leading- and next-to-leading order in 
ChPT. Solid dots and filled rectangles represent vertices from and C^f\ respectively. 



they respect chiral symmetry^] Contrary to dimensionally-regularized expressions, cutoff-regularized 
integrals do not scale properly i.e. their mass dimension is not generated exclusively by the soft scales. 
The renormalized expressions emerging after absorbing the positive powers and logarithms of the cutoff 
into an appropriate redefinition of the LECs do, however, feature the expected scaling behavior. 

I am now in the position to discuss the chiral power counting, i.e. the expression that determines the 
power v of the expansion parameter Q/A x for a given Feynman diagram. This can be achieved by 
carefully counting the powers of small momenta associated with derivatives entering the vertices in C n , 
pion propagators, integrations over the loop momenta and the ^-functions. Using certain topological 
identities, one obtains the following expression for the chiral dimension of a connected Feynman diagram: 

u = 2 + 2L + J2 v i A i, \ = d t -2, (3.41) 

i 

where L refers to the number of loops. This result has been first obtained by Weinberg in [TU]. Notice 
that in order for perturbation theory to work, C n must contain no interactions with Aj < since, 
otherwise, adding new vertices would not increase or even lower the chiral dimension v. This feature 
is guaranteed by the spontaneously broken chiral symmetry of QCD which ensures that only non- 
renormalizable interactions with at least two derivatives or powers of the pion mass appear in C^. In 
particular, chiral symmetry forbids the renormalizable derivative- less interaction of the type 7r 4 . 



Consider now pion-pion scattering as an illustrative example. Eq. (3.41) tells us that the leading 
contribution to the scattering amplitude is generated by tree diagrams (L = 0) constructed from the 
lowest-order vertices with Aj = i.e. the ones from see Fig. [4J The amplitude scales as Q 2 . The 

) grapl 

single insertion from ££ J ', see Fig. 4] 
momenta or one power of the quarl 



(2) 

corrections result from one-loop graphs involving all vertices from C\ as well as tree graphs with a 
,{4) — □ They appear at order Q 4 and are suppressed by two powers of 
rk masses compared to the leading-order contribution. It is easy to 
verify that all diagrams in the bottom line of this figure scale, indeed, as Q 4 . For example, for the first 



J When cutoff regularization is used, a special care is required regarding the treatment of non-covariant pieces in the 
pion propagator, see. |28] . 
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diagram, four powers of momenta arise from the vertices and another four powers of momenta emerge 
from the loop integration. One should further take into account four powers of momenta generated in 
the denominator by the pion propagators. Thus, the total power of the soft scale is indeed four. All 
ultraviolet divergences entering the loop integrals are local and absorbable into redefinition of the LECs 
in cffl (when using dimensional regularization) , as it represents the most general, approximately chiral 
invariant, local interaction of Goldstone bosons at order Q 4 . The divergent parts of the LECs U have 
been worked out in [20] using the heat-kernel method. The finite parts of the l^s are not fixed by chiral 
symmetry and have to be determined from the data or lattice QCD calculations. In the Goldstone 
boson sector, even a number of two-loop calculations (i.e. at order Q 6 ) have already been performed, 
see Ref. [29J for a review article. In particular, an impressive theoretical prediction has been made 
for the isoscalar S-wave ttit scattering length a[j by Colangelo et al. [30J who combined the two-loop 
calculation |31| with dispersion relations to predict a{] = 0.220 ± 0.005. To compare, the leading-order 
calculation by Weinberg yielded aj] = 0.16 [32J while the next-to-leading value obtained by Gasser 
and Leutwyler is ag = 0.20 [20]. The results of the recent E865 experiment at Brookhaven [33J and 
the NA48/2 experiment at CERN [M] combined with the older measurement by the Geneva-Saclay 
collaboration beautifully confirmed the prediction of the two-loop analysis of Ref. [30j yielding the 
value og = 0.217 ± 0.008 (exp) ± 0.006 (th) [35]. This combined result accounts for isospin breaking 
corrections, see Ref. [35] for more details. 

The procedure outlined above can, in principle, be extended to arbitrarily high orders in the low-energy 
expansion. Clearly, the accuracy of the calculations depends crucially on the value of the hard scale A x 
which sets the (maximal) radius of convergence of the chiral expansion. The p-meson is the first meson of 
the non-Goldstone type and shows up as a resonance in p-wave tttt scattering. Such resonances represent 
truly non-perturbative phenomena that cannot be described in standard ChPT^J Consequently, their 
appearance signals the breakdown of the chiral expansion. This leads to the estimation A x ~ M p ~ 770 
MeV. A related observation that matches naturally the above estimation was made by Manohar and 
Georgi who pointed out that A x cannot be larger than AitF^ ~ 1200 MeV since this number sets the 
scale that controls the running of the renormalized LECs when shifting the renormalization point. 

Last but not least, I would like to summarize and underline the special role and importance of the 
chiral symmetry for the whole approach. First of all, it implies severe constraints on the interactions 
in the effective Lagrangian and relates the strengths of various multi-pion vertices. For example, the 



(2) 

leading-order Lagrangian C\ in Eq. (3.36) gives rise to infinitely many vertices, when expanded in 



(2) 

powers of pion fields, whose strengths are determined by just two (!) LECs F and B. CX allows 
to compute the leading contribution to scattering amplitudes for multi-pion processes and to relate 
the strengths of the corresponding matrix elements, thus featuring a remarkable predictive power. 
Moreover, as we saw through the explicit construction, the spontaneously broken chiral symmetry of 
QCD prevents the appearance of derivative-less interactions between pions in the effective Lagrangian. 
The only derivative- less interactions in C n are due to explicit chiral symmetry breaking in Cqcd 
and are suppressed by powers of the light quark masses. When calculating S-matrix elements, the 
derivatives entering the vertices generate powers of external momenta. Consequently, the interaction 

6 For an extension of ChPT to the resonance region, the so-called unitarized ChPT, see |36| and references therein. 
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between pions becomes weak at vanishingly low energies and would even completely disappear if chiral 
symmetry were exact. This turns out to be a general feature of Goldstone bosons and is not restricted 
to the SU (2) l x SU (2)# group. This allows to compute low-energy hadronic observables in a systematic 
way via the chiral expansion, i.e. the dual expansion in powers of momenta and quark masses about 
the kinematical point corresponding to the free theory (assuming that the actual quark masses in the 
real world are low enough for such an expansion to converge). 



3.4 Inclusion of nucleons 

So far we only discussed interactions between Goldstone bosons. We now extend these considerations 
to include nucleons. More precisely, we are interested in describing reactions involving pions with 
external momenta of the order of M n and (essentially) non-relativistic nucleons whose three-momenta 
are of the order of M n . Similarly to the triplet of pion fields, the isospin doublet of the nucleon fields 
should transform nonlinearly under the chiral SU(2)l x SU(2)r but linearly under the vector subgroup 



SU (2)y. The unitary matrix U introduced in Eq. (3.26) is less useful when constructing the Lagrangian 
involving the nucleons. It is more convenient to introduce its square root u, U = u 2 . The transformation 



properties of u under chiral rotations can be read off from Eq. (3.27): 



u — y u = VLuItf = LubT 1 = huR j , (3.42) 
where I have introduced the unitary matrix h = h(L, R, U) given by h = VWrJ LVU which 



is 



sometimes referred to as a compensator field. The last equality in Eq. (3.42) follows from U' = u'u' 
Luh~^u' = LuuR^ . Notice that since pions transform linearly under isospin rotations corresponding to 
L = R = V with U -> U' = VUV j and, accordingly, u — > v! = VuV\ the compensator field in this 
particular case becomes [/-independent and coincides with V. 

Exercise: calculate the explicit form of the compensator field h(L, R, tv) for infinitesimal chiral 



transformations using Eq. (3.32) and keeping only terms that are at most linear in the pion fields. 



Verify that h indeed reduces to the isospin transformation for L = R = V . 

It can be shown that {U, N} define a nonlinear realization of the chiral group if one demands that 

N -»• N' = hN . (3.43) 

I do not give here the proof of this statement and refer the interesting reader to Ref. |23|, 124] . Moreover, 
this nonlinear realization obviously fulfills the desired feature that pions and nucleons transform linearly 
under isospin rotations. Similarly to the purely Goldstone boson case, one can show that all other 
possibilities to introduce the nucleon fields are identical with the above realization modulo nonlinear field 
redefinitions. The most general chiral invariant Lagrangian for pions and nucleons can be constructed 
from covariantly transforming building blocks, i.e. 0% — > 0\ = hOih~ l , by writing down all possible 
terms of the form NO\ . . . O n N. The covariant (first) derivative of the pion field is given by 

U>1 = i«t(0 M i7) u t = -1-^ + 0(tt 3 ) vf = hu^h- 1 , (3.44) 
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Figure 5: Leading pion loop contribution to the nucleon self energy. Solid line represents the nucleon. 



and is sometimes referred to as chiral vielbein. The derivative of the nucleon field, d^N, does not 
transform covariantly, i.e. d^N — > (d^N)' ^ hd^N since the compensator field h does, in general, 
depend on space-time (through its dependence on U). The covariant derivative of the nucleon field 
D„N, D^N -> (iyV)' = hD„N, is given by 

D„N = (d„ + TJN , with = - (Vc> + ud^ = ^ r • tt x d^Tv + 0(tt 4 ) . (3.45) 

The so-called connection can be used to construct higher covariant derivatives of the pion field, for 
example: 

Ufiu = df^Uy + [T M , u„] . (3.46) 
To first order in the derivatives, the most general pion-nucleon Lagrangian takes the form |37j 

4'jv = N (V ^ " m + Y7"75^) iV , (3.47) 

where m and are the bare nucleon mass and the axial-vector coupling constant and the superscript 
of C^N denotes the power of the soft scale Q. Contrary to the pion mass, the nucleon mass does 
not vanish in the chiral limit and introduces an additional hard scale in the problem. Consequently, 



terms proportional to Dq and m in Eq. (3.47) are individually large. It can, however, be shown that 
(i'y^D^ — m)N ~ 0(Q) [3H]. The appearance of the additional hard scale associated with the nucleon 
mass invalidates the power counting for dimensionally regularized expressions since the contributions 
from loop integrals involving nucleon propagators are not automatically suppressed. To see this consider 
the correction to the nucleon mass mjv due to the pion loop shown in Fig. [5] Assuming that the nucleon 
and pion propagators scale as 1/Q and 1/Q 2 , respectively, and taking into account Q 4 from the loop 
integration and Q 2 from the derivatives entering the ^-vertices, the pion loop contribution to the 
nucleon self energy is expected to be of the order ~ Q 3 . Consequently, the corresponding nucleon 
mass shift 5m^ = £(m7v) is expected to be oc M% (since no other soft scale is left). Explicit calculation, 
however, shows that the resulting nucleon mass shift does not vanish in the chiral limit |37j : 

^| loop , rel - 3 §P 3 + ^ In ^) + 0(d - 4) , (3.48) 



where the quantity L(/j,) is defined in Eq. (3.40). The result in Eq. ( |3.48 ) implies that the nucleon mass 



receives a contribution which is formally of the order ~ m (m/47ri ? p~and is not suppressed compared 
to m. The bare nucleon mass m that enters the lowest-order Lagrangian C^L gets renormalized. This 
is in contrast to the purely mesonic sector where loop contributions are always suppressed by powers of 
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(2) 

the soft scale and the parameters F and B in the lowest-order Lagrangian L\ remain unchanged by 
higher-order corrections (if mass-independent regularization is used). I emphasize, however, that even 
though DR expressions do not automatically obey the dimensional power counting with nucleons being 
treated relativistically, the proper scaling in agreement with naive dimensional analysis can be restored 
via appropriately chosen renormalization conditions [39 . Stated differently, one can (and should in 
order for the EFT to be useful) choose renormalization conditions in such a way, that all momenta 
flowing through diagrams are effectively of the order of Q. Another, simpler way to ensure the proper 
power counting exploits the so-called heavy-baryon formalism [4U(, |4"T] which is closely related to the 
nonrelativistic expansion due to Foldy and Wouthuysen [42] and is also widely used in heavy-quark 
effective field theories. The idea is to decompose the nucleon four- momentum according to 

P M = mvp + k^ , (3.49) 

with the four- velocity of the nucleon satisfying v 2 = 1 and fe« its small residual momentum, v-k <C m. 
One can thus decompose the nucleon field N in to the velocity eigenstates 

N v = e imv - x P+N , h v = e imv - x p-N , (3.50) 

where = (1 ± 7 At u M )/2 denote the corresponding projection operators. In the nucleon rest-frame 
with Vfj, = (1,0,0,0), the quantities N v and h v coincide with the familiar large and small components 
of the free positive-energy Dirac field (modulo the modified time dependence). One, therefore, usually 
refers to N v and h v as to the large and small components of N . The relativistic Lagrangian in 



Eq. (3.47) can be expressed in terms of N v and h v as: 

4ii = NvAN v + h v BN v + N vlo B j Joh v - h v Ch v , (3.51) 

where 

A = i(v • D) + g A (S ■ u) , B = -7 5 2i(S ■ D) + ~{v ■ u) , C = 2m + i(v • D) + g A (S • u) , (3.52) 

and = i~f5cr flu v L ' is the nucleon spin operator. One can now use the equations of motion for the 
large and small component fields to completely eliminate h v from the Lagrangian. Utilizing the more 
elegant path integral formulation [43] . the heavy degrees of freedom can be integrated out performing 
the Gaussian integration over the (appropriately shifted) variables h v , h v . This leads to the effective 
Lagrangian of the form [4T] 



A + ( 7 oi3 t 7o)C- 1 el N v = N v [i(v ■ D) + g A (S ■ u)} N v + O ( -) . (3.53) 
J \m J 

Notice that the (large) nucleon mass term has disappeared from the Lagrangian, and the dependence on 
m in C C ^ N resides entirely in new vertices suppressed by powers of 1/m. The heavy-baryon propagator of 
the nucleon is simply 1 / (v ■ k + ie) and can be obtained from the 1/m expansion of the Dirac propagator 



using Eq. (3.49) and assuming v ■ k <C m: 



*fm = A±^ +0(m -, } (3 . 54) 
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where A + = (j) + m) / (2m) is a projection operator on the states of positive energy. The advantage of 
the heavy-baryon formulation (HBChPT) compared to the relativistic one can be illustrated using the 
previous example of the leading one-loop correction to the nucleon mass 



<5mAr| 



loop, HB 



32vrF 2 



(3.55) 



Contrary to the relativistic CHPT result in Eq. (3.48), the loop correction in HBChPT is finite (in DR) 



and vanishes in the chiral limit. The parameters in the lowest-order Lagrangian do not get renormalized 



due to higher-order corrections which are suppressed by powers of Q/A x . Notice further that Eq. (3.55) 
represents the leading contribution to the nucleon mass which is nonanalytic in quark masses. It agrees 



with the result obtained by Gasser et al. based on the relativistic Lagrangian in Eq. (3.47) |37j . In 



general, the power v of a soft scale Q for connected contributions to the scattering amplitude can be 



read off from the extension of Eq. (3.41) to the single-nucleon sector which has the form: 



with 



A.; 



-2 + \m { + di , 



(3.56) 



with rii being the number of nucleon field operators at a vertex i with the chiral dimension Aj. Notice 
that no closed fermion loops appear in the heavy-baryon approach, so that exactly one nucleon line 
connecting the initial and final states runs through all diagrams in the single-baryon sector. 

The heavy-baryon formulation outlined above can be straightforwardly extended to higher orders in 
the chiral expansion. At lowest orders in the derivative expansion, the effective Lagrangian £ Ai for 
pions and nucleons takes the form 
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E 



(NN){NtN) ■ (NtN) 



(3.57) 



Here, the ellipses refer to terms which do not contribute to the nuclear forces up to next-to-next-to- 
leading order (N 2 LO) except for where I have shown only a few terms in order to keep the 
presentation compact. Further, here and in what follows I omit the subscript v of the nucleon field 
operators. The quantity x+ = u^x u ^ + ux^u with x = 2BA4 involves the explicit chiral symmetry 
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breaking due to the finite light quark masses and = O — (0}/2. Finally, Cj, d{, Ci, Ci, D and E 
denote the corresponding LECs. 

The presented elementary introduction into ChPT aims at providing the main conceptual ideas of this 
framework and is neither complete nor comprehensive. Excellent lecture notes on the discussed and 
related subjects |45H5U| are highly recommended for further reading. A very comprehensive, textbook- 
like lecture notes can be found in Ref. [51J. Current frontiers and challenges in these fields are addressed 
in recent review articles [HH [53] , see also Ref. [5l] . 



4 EFT for two nucleons 

4.1 ChPT and nucleon-nucleon scattering 

As outlined in the previous section, ChPT can be straightforwardly extended to the single-nucleon 
sector (apart from the complication related to the treatment of the nucleon mass). A generalization to 
processes involving two and more nucleons is much more difficult. Contrary to the interaction between 
Goldstone bosons, nucleons do interact with each other in the limit of vanishingly small momenta and 
quark masses. Chiral symmetry does not constrain few-nucleon operators in the effective Lagrangian 



which contains derivative-less terms, see Eq. (3.57). In fact, the interaction between the nucleons at 



low energy is even strong enough to bind them together. Shallow bound states such as the deuteron, 
triton etc. represent non- perturb at ive phenomena that cannot be described in perturbation theory. 

On the other hand, just following the naive dimensional analyses as we did in the previous section, 
the power counting can be straightforwardly generalized to connected Feynman diagrams involving N 
nucleons leading to 

v = 2- N + 2L + Y^ ViAi . (4.58) 

i 

This implies the usual suppression for loop diagrams and thus suggests that the interaction is weak. 
This conclusion is certainly not correct. So, what goes wrong? The reason why the naive dimensional 
analysis yields a wrong result is due to the appearance of infrared divergences (in the HBChPT) 
in diagrams which contain purely nucleonic intermediate states [U [5]. Consider the two-pion- (2ir-) 
exchange box Feynman diagram shown in Fig. [6] (the diagram on the left-hand side). In the nucleon 
rest frame with u„ = (1, 0, 0, 0), the four-momenta of the incoming nucleons are (k 2 / (2m) + 0(m~ 3 ), k) 
and (k 2 /(2m) + 0(m~ 3 ), —k). In the infrared regime with k = 0, the contribution of the box diagram 
takes the form 

dH P(l) 



(2vr) 4 (JO + ie)(-l° + ie)(l 2 - M 2 + ie) 2 ' (4 " 59) 

where I is the loop momentum and P(l) is a polynomial whose explicit form is determined by the 
pion-nucleon vertex. The integral over 1° possesses the so-called pinch singularity due to the poles at 
1° = ±ie. Notice that such pinch singularities only show up in the case of at least two nucleons since for 
a single nucleon the contour of integration can be distorted to avoid the singularity. The singularities 
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(a) (b) (c) (d) (e) (f) (g) 

Figure 6: Two-pion exchange: Feynman diagram (a) and the corresponding time-ordered graphs (b-g). 
Solid (dashed) lines correspond to nucleons (pions). 



that appear in the box diagram are not "real" but an artefact of the heavy-baryon approximation for 
the nucleon propagators (static nucleons) that is not valid for such diagrams. 

Exercise: verify this statement by using the standard Dirac propagators for the nucleon field, see 
Eq. (3.54), an d making the nonrelativistic expansion after carrying out the integration over 1°. 



An alternative and, perhaps, more instructive way to explore the origin of the infrared enhancement is 
by using the so-called "old-fashioned" time-ordered perturbation theory instead of the covariant one. 
In time-ordered perturbation theory, the T-matrix is given by 



tu\ , ST^ ( H l)aa(Hl) a i3 ^ (#/)aa(#/)afe(-ft/)&/3 ( a cn\ 

^ Ep - E a + ie (Ep - E a + ie){Ep - E b + it) 



where Hj is the interaction Hamiltonian corresponding to the effective Lagrangian for pions and nu- 
cleons. This expression should be familiar from Quantum Mechanics. Its derivation and application to 
quantum field theory can be found e.g. in [55] ■ Here, I use Latin letters for intermediate states, which, 
in general, may contain any number of pions, in order to distinguish them from purely nucleonic states 
denoted by Greek letters. I remind the reader that no nucleon-antinucleon pairs can be created or 
destroyed if nucleons are treated nonrelativistically. Consequently, all states contain the same number 
of nucleons. It is useful to represent various contributions to the scattering amplitude in terms of 
time-ordered diagrams. For example, the Feynman box diagram for NN scattering via 27r-exchange can 
be expressed as a sum of six time-ordered graphs, see Fig. [6j which correspond to the following term 



in Eq. (4.60): 



E (H 7T NN)aa(H 7T NN)ab(H 7T NN)bc(H 7T NN)cl3 (A (11, 

ahc (Ep -E a + ie)(Ep -E b + ie)(Ep -E c + ie) ' 1 ' ' 

where H n jyj^ denotes the ttNN vertex. Actually, this expression can be obtained from carrying out the 
/°-integration in the corresponding Feynman diagram (using Dirac propagators for the nucleons). It is 
easy to see that the contributions of diagrams (d-g) are enhanced due to the presence of the small (of 
the order Q 2 /m) energy denominator associated with the purely nucleonic intermediate state \b) which 
in the center-of-mass system (CMS) takes the form: 

(4.62) 



Ep - E b + ie Pp/m - p^/m + ie 
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Notice that the energy denominators corresponding to the ttNN states \a) and |c) contain the pion 

energy uj^ = \Jk 2 + M% and are of the order M w ~ Q in agreement with the dimensional analysis. 
According to Weinberg, the failure of perturbation theory in the few-nucleon sector is caused by the 
enhanced contributions of reducible diagrams, i.e. those ones which contain purely nucleonic intermedi- 
ate states. It should, however, be emphasized that the infrared enhancement is not sufficient to justify 
the need of non-perturbative resummation of the amplitude if one counts m = 0(A X ). According 
to Eq. (4.58) and taking into account the infrared enhancement ~ m/Q due to the purely nucleonic 
intermediate states, loop contributions are still suppressed by ~ Qm/A^ ~ Q/A for m ~ A x . To 
overcome this conceptual difficulty, Weinberg proposed to treat the nucleon mass as a separate hard 
scale according to the rule [U [5] : 

A? 



m ~ 



Q 



> A. 



(4.63) 



The resulting power counting is referred to as the Weinberg power counting. I will also discuss some 
alternative scenarios. 

The infrared enhancement of the few-nucleon diagrams can be naturally taken into account by re- 



arranging the expansion in Eq. (4.60) and casting it into the form of the Lippmann-Schwinger (LS) 
equation 



'-a/3 



effJa/3 



E, 



E~ + ie ' 



(4.64) 



with the effective potential (V e s) a i3 defined as a sum of all possible irreducible diagrams (i.e. the ones 
which do not contain purely nucleonic intermediate states): 



eS)a/3 



(Hl)aa(Hl)aj3 

E B -E a + ie 



(Hl)aa(Hl)ab(H, 



I)bj3 



ab 



{Efs-E a + ie){E p -E b + ie] 



+ 



(4.65) 



Here, the states |a), \b) contain at least one pion. The effective potential in Eq. (4.65) does not contain 



small energy denominators and can be worked out within the low-momentum expansion following the 
usual procedure of ChPT. After the potential is obtained at a given order in the chiral expansion, 



few-nucleon observables can be computed by solving the LS equation (4.64), which leads to a nonper 



turbative resummation of the contributions resulting from reducible diagrams. The resulting two-step 
approach will be referred to as ChEFT in order to distinguish it from ChPT in the Goldstone boson 
and single-nucleon sectors. 



4.2 Analytic properties of the non-relativistic scattering amplitude 

Before discussing various scenarios of organizing EFT for two nucleons, it is useful to recall general con- 
straints imposed on the partial wave scattering amplitude by analyticity. Consider two non-relativistic 
nucleons interacting via a potential V. The corresponding S-matrix for an uncoupled channel with the 
orbital angular momentum / is parametrized in terms of a single phase shift 5i and can be written in 
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Figure 7: Some singularities of the partial wave T-matrix (left panel), effective range function Fi(E) 
(middle panel) and the modified effective range function F l M (E) (right panel) . The shaded areas show 
the (maximal) range of applicability of the ERE and MERE. 



terms of T-matrix as 

S l = e^ = l- i (^T l (k), (4.66) 

with k denoting the CMS scattering momentum. The T-matrix can then be expressed in terms of the 
so-called effective range function = k 2l+l cot5i(k) via 

16^ 2 k 21 

w = -VfipH^- (467) 

In the complex energy plane, the scattering amplitude and thus also the T-matrix possess a so-called 
unitarity cut, a kinematic singularity due to two-body unitarity. The unitarity cut starts from the 
branch point at the threshold (E = 0) and goes to positive infinity. The dynamic singularities are 
associated with the interaction mechanism and are located at the negative real axis. For example, in 
the case of Yukawa potential ~ exp(— Mr)/r corresponding to an exchange of a meson of mass M, the 
amplitude has a left-hand cut starting at k 2 = —M 2 /4. Bound and virtual states reside as poles at 
the negative real axis (k = i\k\ and k = —i\k\ for bound- and virtual-state poles, respectively) while 
resonances show up as poles at complex energies. 

Exercise: verify the appearance of the left-hand cut in the scattering amplitude for the one-pion exchange 
potential 

V(j)', p) oc °\ , with q = p'-p, 

using the first Born approximation. 

Contrary to the scattering amplitude, the effective range function does not possess the kinematic 
unitarity cut and can be shown to be a real meromorphic function of k 2 near the origin for non-singular 
potentials of a finite range [561 E7|- It can > therefore, be Taylor-expanded about the origin leading to 
the well-known effective range expansion (ERE) 

F^k 2 ) = -- + Kk 2 + v 2 k 4 + v 3 k 6 + . . . , (4.68) 
a 2 
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with a, r and Vj being the scattering length, effective range and the so-called shape parameters. Gen- 
erally, the maximal radius of convergence of the ERE is limited by the lowest-lying left-hand dynamic 
singularity associated with the potential. For Yukawa-type potentials with the range r ~ M~ x , the 
(maximal) radius of convergence of the ERE is given by k 2 < M 2 /4. For (strong) nucleon-nucleon 
interaction with the one-pion- (lvr-) exchange potential constituting the longest-range contribution, the 
ERE is expected to converge for energies up to ~ M 2 / (2rriN) = 10.5 MeV. Notice that apart from 
the singularities associated with the structure of the potential, Fi{k 2 ) may also contain discrete poles 
whose positions are determined by the strength of the interaction. The appearance of such poles near 
the origin would spoil the convergence of the ERE. 

The framework of ERE can be generalized to the case in which the potential is given by a sum of a long- 
range (r/ ~ m^ 1 ) an d short-range (r s ~ m^ 1 <C mj ) potentials Vl and Vg, respectively. Following 
van Haeringen and Kok [58], one can define the modified effective range function F, via 



F t M (k 2 ) = Mt(k) + jjj^ cotft(*) - St(k)\ , (4.69) 



The Jost function //'(/e) is defined according to f("{k) = fi L {k,r)\ r=Q with f["{k,r) being the Jost 
solution of the Schrodinger equation corresponding to the potential Vl, i-e. the particular solution that 
fulfills 

lim e~ ikr fi{k, r) = 1. (4.70) 



r— >co 



Further, 5f{k) denotes to the phase shift associated with the potential Vl and the the quantity M^{k) 
can be computed from f^{k,r) as follows: 



MHk) = lim 



I <t~ 1 ; III /%' /' ! 

(4.71) 



ik\' i .. r d 2l+1 ,/,(/,•. r 



r 



dr 2l+l f^ty 



I denote here with the superscript "L" all quantities that can be computed solely from the long-range 
part of the potential. The modified effective range function F^ (k 2 ) defined in this way does not contain 
the left-hand singularity associated with the long-range potential and reduces, per construction, to the 
ordinary effective range function Fi(k 2 ) for V L = 0. It is a real meromorphic function in a much larger 
region given by rj 1 as compared to i^(/c 2 )(^] If the long-range interaction is due to a Coulomb potential, 
Vl{t) = a/r, the Jost solution and, consequently, the function M^(k) can be calculated analytically. 
For example, for I = and the repulsive Coulomb potential, the MERE takes the following well-known 
form: 

F c {k 2 ) = Cl{r,) k cot[8{k) - 5 C (k)} + 2& 77/1(77) , (4.72) 
where the Coulomb phase shift is 6 = arg r(l + irj) and the quantity r\ is given by 

Tfi 

77= -a. (4.73) 



7 Note that the existence of Mj \(k) implies certain constraints on the small-r behavior of Vz,(r). 
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Further, the functions Cq(t)) (the Sommerfeld factor) and h(rj) read 

-12 ,^ 27T7 1 



o2irri _ I 



and h{rf) = Re ty(irj) — hi{rf) 



(4.74) 



Here, ^f(z) = T'(z)/T(z) denotes the digamma function. For more details on the analytic properties of 
the scattering amplitude and related topics I refer the reader to the review article |59j . 

After these preparations, we are now in the position to discuss the implications of the long-range 
interaction on the energy dependence of the phase shift. It is natural to assume that the coefficients in 
the ERE and MERE (except for the scattering length) are driven by the scales mi and m s associated 
with the lowest left-hand singularities, see [60J for a related discussion. The knowledge of the long-range 
interaction Vl allows to compute the quantities fi L (k), M t L (k) and 5f{k) entering the right-hand side 
of Eq. (4.69) and thus to express Si(k) and the ordinary effective range function Fi(k 2 ) in terms of the 
modified one, F t M (k 2 ). The MERE for F l M (k 2 ) then yields an expansion of the subthreshold parameters 
entering Eq. (4.68) in powers of mi/m s . In particular, using the first few terms in the MERE as input 



allows to make predictions for all coefficients in the ERE. The appearance of the correlations between 
the subthreshold parameters in the above-mentioned sense which I will refer to as low-energy theorems 
(LETs) is the only signature of the long-range interaction at low energy (in the two-nucleon system). 
The LETs allow to test whether the long-range interactions are incorporated properly in nuclear chiral 
EFT and thus provide an important consistency check. 



4.3 EFT for two nucleons at very low energy 

Before discussing chiral EFT for two nucleons, let us consider, as a warm-up exercise, a simpler EFT 
for very low energies with Q <C M n . Then, no pions need to be taken into account explicitly, and the 
only relevant degrees of freedom are the nucleons themselves. The corresponding EFT with the hard 
scale A ~ M n is usually referred to as pionless EFT. The most general effective Lagrangian consistent 
with Galilean invariance, baryon number conservation and the isospin symmetry takes in the absence 
of external sources the following form: 

/ V 2 \ 1 1 

C = NUid + ) N — -C s (At AT) (At AT) - -C T (N^aN) ■ (A^A) + . . . , (4.75) 

\ 2m 12 2 

where Cst are LECs and the ellipses denote operators with derivatives. Isospin-breaking and relativistic 



corrections to Eq. (4.75) can be included perturbatively. 

What can be expected from the pionless EFT as compared to the ERE? In the absence of external 
sources and restricting ourselves to the two-nucleon system, both approaches provide an expansion of 
NN low-energy observables in powers of k/M n , have the same validity range and incorporate the same 
physical principles. Pionless EFT can, therefore, not be expected to do any better than ERE. Our goal 
will be thus to design the EFT in such a way that it matches the ERE for the scattering amplitude 

16tt 2 1 

(4.76) 



m 



■\ + \r k 2 + v 2 k A + v 3 k 6 + ...)- ik 
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Here, I restrict myself to S-waves only. While the coefficients in the effective range expansion are, in 
general, driven by the range of the potential and thus expected to scale with the appropriate powers of 
M^, the scattering length can, in principle, take any value. In particular, it diverges in the presence of 
a bound or virtual state at threshold. It is, therefore, useful to distinguish between a natural case with 
\a\ ~ M" 1 and an unnatural case with \a\ ^> M" 1 . In the natural case, the T-matrix in Eq. (4.76) can 
be expanded in powers of k as: 



16n 2 a 



in 



1 



iak + ( — - 
V 2 



k 2 + 



(4.77) 



where the superscripts of T denote the power of the soft scale Q. A natural value of the scattering 
length implies the absence of bound and virtual states close to threshold. The T-matrix can then be 
evaluated perturbatively in the EFT provided one uses an appropriate renormalization scheme (i.e. the 
one that does not introduce an additional large scale). When solving the LS equation with point-like 
contact interactions, one encounters divergent loop integrals of the kind 
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with n = 1, 3, 5, . . . , I(k) 
The integrals can be evaluated using a cutoff regularization: 
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(4.78) 
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(4.79) 



where the last equation is valid for k < A. To renormalize the scattering amplitude, I divide loop 
integrals into the divergent and finite parts and take the limit A — > od2 



lim I, 

A— »oo 



A - lim f/A + !!^ 
A^oo \ n 2mr 2 



2nvr 2 
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nifj, 
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imk 



A(f,)+I R (ti,k) 



(4.80) 



Here, A n (fi n ) and A(/x) denote the divergent parts of the loop integrals while I„ (/i n ) and I R (fi,k) are 
finite. The procedure is analogous to the standard treatment of divergences arising from pion loops 
in ChPT, see section [5J The splitting of loop integrals in Eq. (4.80) is not unique. The freedom in 
the choice of renormalization conditions is parameterized by \x and fx n . The divergent parts A n (/x n ) 
and A(/x) are to be canceled by contributions of counterterms. The renormalized expression for the 



amplitude, therefore, emerge from dropping the divergent parts in Eq. (4.80) and replacing the bare 
LECs by the renormalized ones Cj — > C^({/i, /U n }). The proper choice of renormalization conditions 



8 While extremely convenient in the case under consideration, taking the limit A — > oo is, strictly speaking, not necessary 
in an EFT. It is sufficient to ensure that the error from keeping the cutoff finite is within the accuracy of the EFT expansion. 
In the considered case, taking A ~ would do equally good job in describing the scattering amplitude. 
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Figure 8: Leading, subleading and sub-subleading contributions to the S-wave T-matrix in the case of 
a natural scattering length. Solid dots (filled rectangles) refer to contact vertices without (with two) 
derivatives. Lines represent the nucleon propagators. 



requires choosing fx, fi n ~ Q <C M n . Dimensional regularization with the minimal subtraction can be 
viewed as a special case corresponding to fi = fj, n = 0. Another special case is given by DR with the 
power divergence subtraction (PDS) [61, 62J. In this scheme, the power law divergences, which are 
normally discarded in DR, are explicitly accounted for by subtracting from dimensionally regulated 
loop integrals not only l/(d— 4)-poles but also the l/(d— 3)-poles. Its formulation used in Refs. [HUES] 
corresponds to the choice fi n = and \i — > fiir/2 in Eq. ( 4.80| ). 

The dimensional analysis for the renormalized scattering amplitude implies that the leading and sub- 
leading terms and T^ 1 ) are given by the tree- and one-loop graphs constructed with the lowest-order 
vertices from Eq. (4.75), see Fig.g T( 2 ) receives a contribution from the two-loop graph with the lowest- 
order vertices and the tree graph with a subleading vertex [61\ 162] . Higher-order corrections can be 



evaluated straightforwardly. Matching the resulting T-matrix to the ERE in Eq. (4.77) order by order 
in the low- momentum expansion allows to fix the LECs C\. At next-to-next-to- leading order (N 2 LO), 
for example, one finds: 



m l 



1 + O(afi) 



2-rra 2 



r . 



(4.81) 



where the LECs Co and C2 are defined vict th.6 tree-level X'-nicitrix: 2~tree 
LEC C is related to C S t in Eq. (14. 751) as C = C s 



4vr(C + C 2 fc 2 + ...). The 
3Cr- Here and in the remaining part of this 
section, the expressions are given in DR with PDS. 

For the physically interesting case of neutron-proton scattering, the two S-wave scattering lengths 
appear to be large: 



r/i 



So 



-23.714 fm 



-16.6 M. 



-1 



r/.-i 



Si 



5.42 fm~ 3.8 M: 



(4.82) 



Instead of using the low-momentum representation in Eq. (4.77) which is valid only for k < 1/a, it is 
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'(-1) = 



, (0) 



where: 



Figure 9: Leading and subleading contributions to the S-wave T-matrix in the case of an unnaturally 
large scattering length. For notation, see Fig. [8} 



advantageous to expand the T-matrix in powers of k keeping ak ~ 1 [6T| 152]: 
T = T^ 1 ) + T^ ) + T^ 1 ) + . . . 



(4.83) 
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The EFT expansion can be cast into the form of Eq. (4.83) if one assumes that the LECs CT scale 
according to C2n ~ 1/Q n+1 . The leading contribution T 1 " 1 * 1 then results from summing up an infinite 
chain of bubble diagrams with the lowest-order vertices, see Fig. [9} All diagrams constructed only from 
Cq scale as 1/Q. For example, for the one- loop graph one has Q 3 from the integration, 1/Q 2 from the 
nucleon propagator and another 1/Q 2 from the LECs Cq. The corrections are given by perturbative 
insertions of higher-order interactions dressed to all orders by the leading vertices. Matching the 



resulting T-matrix with the one in Eq. (4.83) one finds at NLO: 
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(4.84) 



Exercise: calculate the S-wave scattering amplitude up to NLO for the case of unnaturally large 



scattering length and verify the expressions for the LECs given in Eq. (4.84)- 



The real power of pionless EFT comes into play when one goes beyond the two-nucleon system by 
considering e.g. low-energy reactions involving external electroweak sources and/or three- and more 
nucleon systems. A discussion of these topics goes beyond the scope of these lectures. I refer an 
interested reader to the recent review articles 1 63 - 66 1 . 
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4.4 Chiral EFT for two nucleons with perturbative pions 



We have seen in the previous section how the EFT without explicit pions can be organized to describe 
strongly interacting nucleons at low energy. The limitation in energy of this approach, cf. the discussion 



in section 4.2 appears to be too strong for most nuclear physics applications. To go to higher energies 
it is necessary to include pions as explicit degrees of freedom. I have already outlined in section [ITT] one 
possible way to extend ChPT to the few-nucleon sector following Weinberg's original proposal [H [5]. 
In this approach, the nonperturbative dynamics is generated by iterating the lowest-order two-nucleon 
potential vS*) which subsumes irreducible (i.e. non-iterative) contributions from tree diagrams with the 



leading vertices (i.e. Aj = 0), see Eq. (4.58). The only possible contributions are due to derivative- less 



contact interaction and the static l7r-exchange, so that the resulting potential reads: 

2N = ~4Ff g2 + M 2 T i -T2 + C s + C T (Jx ■ a 2 . (4.85) 

Here, 5i (r^) are the Pauli spin (isospin) matrices of the nucleon i, q = p' — p is the nucleon momentum 
tran sfer and p (p ') refers to initial (final) nucleon momenta in the CMS. As pointed out in section 
the justification for resumming to all orders in the LS equation is achieved in the Weinberg 



4.1 



approach via a fine tuning of the nucleon mass, see Eq. (4.63). With this counting rule, it follows 



immediately that every iteration of VjjjJ in Eq. (4.64) generates a contribution of the order Q°/A x . On 



the other hand, in the pionless EFT with unnaturally large scattering length outlined in section [4T3| the 
non-perturbative resummation of the amplitude was enforced by fine-tuning the LECs accompanying 
the contact interactions while treating the nucleon mass on the same footing as all other hard scales. 
While these two scenarios are essentially equivalent in the pionless case, they lead to an important 
difference in organizing EFT with explicit pions. The approach due to Kaplan, Savage and Wise 
(KSW) |61] [62] represents a straightforward generalization of the pionless EFT to perturbatively include 
diagrams with exchange of one or more pions. The scaling of the contact interactions is assumed to be 
the same as in pionless EFT (provided one uses DR with PDS or an equivalent scheme to regularize 
divergent loop integrals). Notice that in contrast to pionless EFT, the pion mass is treated as a soft 
scale with Q ~ M n ~ a -1 . The only new ingredients in the calculation of the amplitude up to next-to- 
leading order in the KSW expansion are given by the dressed l7r-exchange potential and derivative-less 



interaction cx M 2 , see Fig. 10 27r-exchange is suppressed and starts to contribute at N 2 LO. At each 
order in the perturbative expansion, the amplitude is made independent on the renormalization scale 
by an appropriate running of the LECsj^] Compact analytic expressions for the scattering amplitude 
represent another nice feature of the KSW approach. 



As explained in section 4.2 the appearance of a long-range interaction implies strong constraints on 
the energy dependence of the amplitude and imposes certain correlations between the coefficients in 
the ERE (LETs). EFT with explicit pions aims at a correct description of non-analyticities in the scat- 
tering amplitude associated with exchange of pions which in this framework represent truly long-range 
phenomena. Thus, the correct treatment of the long-range interaction by including pions perturbatively 



9 Strictly speaking, an exact scale independence of the NLO amplitude in the KSW approach with explicit pions is 
achieved at the cost of resumming a certain class of higher-order terms, see the discussion in Ref. |67j . 
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Figure 10: Leading and subleading contributions to the S-wave T-matrix in the case of unnaturally 
large scattering length in the KSW approach with explicit pions. Filled rectangle denotes contact 
interactions with two derivatives or a single insertion of M 2 . For remaining notation, see Fig. |y 



can be ultimately judged by testing the corresponding low-energy theorems. This idea was picked up by 
Cohen and Hansen [68, 69J. Fixing the LECs accompanying the contact interactions from the scatter- 
ing length and effective range, they obtained the following predictions for the first three S-wave shape 
parameters at NLO in the KSW scheme [68, 69J: 

g\m ( 16 32 2 \ 



V2 16vrF 2 V 3a 2 Af* + 5aM% M 2 J 
g\m ( 16 128 16 



16vrF 2 \a?M% 7aM% 3M* 
g\m ( 256 | 512 16 \ 
VA ~ 16^Ff V 5oJm! + 9aJIJ ~ Ml ) " ( } 

Plugging in the numbers for the nucleon and pion masses, qa — 1-27, = 92.4 MeV and the corre- 
sponding scattering lengths, see Eq. (4.82), Cohen and Hansen obtained the results quoted in Table l| 10 



A clear failure for the LETs in both channels serves as an indication that the long-range physics is not 
properly taken into account if pions are treated perturbatively. The convergence of the KSW expansion 
was further tested in Ref. |71j where the amplitude is calculated up to N 2 LO. While the results for 
the 1 Sq and some other partial waves including spin-singlet channels were found to be in a reasonable 
agreement with the Nijmegen PWA, large corrections show up in spin-triplet channels already at mo- 
menta of the order of ~ 100 MeV and lead to strong disagreements with the data. The perturbative 
inclusion of the pion-exchange contributions does not allow to increase the region of validity of the 
EFT compared to the pionless theory, see however, Ref. {72j for a new formulation which is claimed 
to yield a convergent expansion. The failure of the KSW approach in the spin-triplet channels was 
attributed in pT] to the iteration of the tensor part of the l7r-exchange potential. This appears to be in 
line with phenomenological successes of Weinberg's approach which treats pion exchange contributions 
nonperturbatively. The most advanced analyses of the NN system at next-to-next-to-next-to-leading 
order (N 3 LO) in Weinberg's power counting scheme demonstrate the ability to accurately describe NN 
scattering data up to center-of-mass momenta at least of the order ~ 2M n [8l I73j. 



3 A slightly different values compared to the ones quoted in this table are extracted from the Nijmegen PWA in Ref. [70] 
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l So partial wave 


3 Si partial wave 




V2 [fm 3 ] V3 [fm 5 ] [fm 7 ] 


V2 [fm 3 ] V3 [fm 5 ] V4 [fm 7 ] 


LETs 
Nijmegen PWA 


-3.3 17.8 -108 
-0.48 3.8 -17 


-0.95 4.6 -25 
0.4 0.67 -4 



Table 1: A comparison of the predicted S-wave shape parameters from Ref. [68J with coefficients 
extracted from the Nijmegen PWA. 



4.5 Towards including pions nonperturbatively: playing with toy models 

While the power counting approach due to Weinberg allows for a nonperturbative resummation of the 
l7r-exchange potential, there is a price to pay. Contrary to the KSW approach, the leading NN potential 
is non-renormalizable in the traditional sense, i.e. iterations of the LS equation generate divergent 
terms with structures that are not included in the original potential. Moreover, resummation of the 
l7r-exchange potential in the LS equation can only be carried out numerically. This prevents the use of 
regularization prescriptions such as e.g. DR which avoid the appearance of a hard scale and maintain the 
manifest power counting for regularized loop contributions making renormalization considerably more 
subtle. Most of the available calculations employ various forms of cutoff regularization with the cutoff 
being kept finite. The purpose of this section is to provide an in-depth discussion on regularization and 
renormalization in this context by considering a simple and exactly solvable quantum mechanical model 
for two nucleons interacting via the long- and short-range forces. This may be regarded as a toy model 
for chiral EFT in the two-nucleon sector. In this model, resummation of the long-range interaction can 
be carried out analytically. This allows to employ and compare the subtractive renormalization that 
maintains the manifest power counting and the cutoff formulation of the effective theory. I also explore 
the consequences of taking very large values of the cutoff in this model. The presentation follows closely 
the one of Ref. [67] . 



4.5.1 The model 



Consider two spin-less nucleons interacting via the two-range separable potential 



V(p', p) = «, F t tf) F,(p) + v s F s (p') F 8 (p) , F{p) 



mi 



p 2 + mf 



F s (p) 



mi 



(4.87) 



where the masses mi and m s fulfill the condition mi -C m s . Further, the dimensionless quantities vi 
and v s denote the strengths of the long- and short-range interactions, respectively. The choice of the 
explicit form of F/ iS (p) is motivated by the simplicity of calculations. 



The potential in Eq. (4.87) does not depend on the angle between p and p and, therefore, only gives 



rise to S-wave scattering. The projection onto the S-wave in this case is trivial and simply yields the 
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factor of 4-7T from the integration over the angles. For an interaction of a separable type, the off-shell T- 
matrix and, consequently, also the coefficients in the ERE, see Eqs. (4.67) and (4.68), can be calculated 
analytically by solving the corresponding LS equation 

T(p', p; k) = V(p', p) + 4vr J j^V(p', l) k2 _ p+i T(l, p; k) , (4.88) 

where m is the nucleon mass and k corresponds to the on-shell momentum which is related to the 
two-nucleon center-of-mass energy via Eqms = k 2 /m. Note that we have absorbed the factor 4tt into 
the normalization of the T-matrix which is, therefore, different from the one in Eq. (2.6). In particular, 
the relation between the S- and T-matrices is given by S(p) = 1 — ipmT(p, p; p)/(2n). 



As explained in section 4.2, the coefficients in the ERE generally scale with the mass corresponding 
to the long-range interaction that gives rise to the first left-hand singularities in the T-matrix. Thus, 
in the considered two-range model, the coefficients in the ERE can be expanded in powers of mi/m s 
leading to the "chiral" expansion: 







1 "a 
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1 "a 
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mf 
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mi 






m s 




m 2 s 





where , a^ 1 and ai™ are dimensionless constants whose values are determined by the form of the 
interaction potential. I fine tune the strengths of the long- and short-range interactions in our model 
in such a way that they generate scattering lengths of a natural size. More precisely, I require that 
the scattering length takes the value a = Q-i/mi (a = a s /m s ) with a dimensionless constant \ai\ ~ 1 
(\a s \ ~ 1) when the short-range (long-range) interaction is switched off. This allows to express the 
corresponding strengths vi and v s in terms of a; and a s as follows: 

87TTO?a/ 47rm,Q!= 

VI = 7 2 , 2 o 2\ ' V ' = r-^V 4 - 90 

m [atmj + mfai — 2mj) m [a s — 1 ) 

One then finds the following expressions for the first three terms in the "chiral" expansion of the 
scattering length 

a£> = ai , = (ai - l) 2 a s , = (ai - l) 2 aia 2 s , (4.91) 



and effective range 



a® 



3ai 



2 (q t - 1) (3a; - 4) a s 



(2 ) (a; - 1) (3a; - 4) (5a; - 3) a 2 s + (2 - a;) a| 



a?* = ^ '-^ ,v \ ' s v ^=*-. (4.92) 

af 
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Notice that in the model considered the leading terms in the ?n//m s -expansion of the ERE coefficients 
are completely fixed by the long-range interaction. The scenario realized corresponds to a strong (at 
momenta k ~ m{) long-range interaction which needs to be treated non-perturbatively and a weak 
short-range interaction which can be accounted for in perturbation theory. 



Exercise: calculate the T-matrix by solving the LS equation (4-88) for the potential given in Eq. (4.81). 
Verify the "chiral" expansion for the scattering length and effective range. Work out the first terms in 
the "chiral" expansion of the shape parameters V2 and V3. The results can be found in Ref. I67f. 



4.5.2 KSW-like approach 

At momenta of the order k ~ mi <C m s , the details of the short-range interaction cannot be resolved. 
An EFT description emerges by keeping the long-range interaction and replacing the short-range one 
by a series of contact terms V^hort(p', p) = Co + C2(p 2 +p' 2 ) + ■ ■ • • Iterating such an effective potential 
in the LS equation leads to ultraviolet divergences which need to be regularized and renormalized. 
The renormalization prescription plays an important role in organizing the EFT expansion. I first 
consider the most convenient and elegant formulation by employing the subtractive renormalization 
which respects dimensional power counting at the level of diagrams. In this sense, the considered 
formulation is conceptually similar to the KSW framework with perturbative pions discussed above 
and will be referred to as the KSW-like approach. The available soft and hard scales in the problem 
are given by Q = {k, /x, mi} and A = {m s , m}, respectively. Here fi ~ mi denotes the subtraction 
point. (There is just a single linearly divergent integral at the order in the low-momentum expansion 
I will consider.) The contributions to the amplitude up to N 2 LO in the Cj/A-expansion are visualized 



in Fig. 11 and can be easily verified using naive dimensional analysis. In particular, the leading term 



arises at order Q 1 is generated by the leading term in the Q/A-expansion of the long-range interaction 
Vion g (p',p) = vi Fi(p') Fi(p) (4-93) 



8-Kmfai 



m (ai — 2) (p 2 + mf)(p' 2 + mf) 



aim 2 jP_ pf_ + (D fQ 4 



(ai - 2) m 2 2m 2 2m 2 \ A 4 



which scales as Cj -1 and, therefore, needs to be summed up to an infinite order, see Fig. 11 Notice that 
the natural size of the short-range effects in our model suggests the scaling of the short-range interactions 
in agreement with the naive dimensional analysis, i.e. C2 n ~ Q ■ This leads to the following expression 
for the on-the-energy shell T-matrix: 

= l — , (4.94) 

m(k — imi) 2 [k 2 (ai — 2) + 2ikm\ (a/ — 2) + 2mJ\ 

from which one deduces 

A:cot<5 = - 47r - 1 ^+^= _^ + (3«z-4) fc2+ (a,-2) ^ 
m T( L > ai 2miai 2mfa>i 
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Figure 11: Leading, next-to- leading and next-to-next-to- leading order contributions to the scattering 
amplitude in the KSW-like approach. The solid lines denote nucleons while the dashed ones represent 
an insertion of the lowest-order (i.e. 0(Q~ 1 )) long-range interaction. Solid dots (dotted lines) denote 
an insertion of the lowest-order contact interaction oc Co (subleading order-0((5) contribution to the 
long-range interaction) . 



Not surprisingly, one observes that the leading terms in the expansion of the ERE coefficients in 
Eq. (4.89) are correctly reproduced. 

The first correction at order Q° is given by the leading-order contact interaction dressed with the 



(4.96) 



iterated leading long-range interaction as visualized in Fig. [TTJ One finds 

C (k + imi) 2 [k 2 (a t - 2) + 2m 2 (at - 1)] 2 



r (o) 



(k - imi) 2 [k 2 (m - 2) + 2ikm t (a/ - 2) + 2m 2 ] 2 ' 



Notice that all integrals entering T( ^ and T(°> are finite. The effective range function k cot 5 at NLO 
can be computed via 



k cot 5 



4vr 1 



in 



1 



r (0) 



+ ik. 



(4.97) 



The "chiral" expansion of the coefficients in the ERE results from expanding the right-hand side in 
this equation in powers of k 2 a nd, subsequently, in powers of m/. The LEC Cq can be determined from 
matching to in Eq. (4.91) which yields 



Co 



Aira s 
mm s 



This leads to the following predictions for the effective range: 

3ai 



1 

mi 



4 2 (ai - 1) (3a/ - 4) a s 

- H o m i 

afm s 



(4.98) 



(4.99) 
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One observes that a?' is correctly reproduced at NLO. The same holds true for the first shape param- 
eters, see Ref. [67] for explicit expressions. Moreover, using dimensional analysis it is easy to verify 
that, in fact, for all i must be reproduced correctly at this order. This is a manifestation of the 
LETs discussed in section l4~2l 

At N 2 LO, the linearly divergent integral I(k) occurs, see Eq. (4.79), which is treated according to 
Eq. (4.80). Renormalization is carried out by dropping its divergent part A(//) and replacing the 
bare LEC Co by the renormalized one Cq(//) in the expression for the amplitude. A straightforward 
calculation using MATHEMATICA yields: 



T (i) 



(k + im{) 



32ir 3 k 2 mf (a, - 2) a L 



Air 2 mm1 (k — imi) 2 [kat[ (k + 2im\) — 2 (k + imi) 2 ] 2 
+ (C (//)) 2 m 2 m 2 s [k 2 (a, - 2) + 2m} (a, - 1)] 2 

a t [k 2 (-2fi - iirk) + 2k{iik - 2i//)m z + 2irmf] + 2(2// + ink) (k + imi) 2 
kai (k + 2im{) — 2{k + imi) 2 

The LEC Co(/t) can be written in terms of a perturbative expansion as follows 

C (m) = cS 0) + C«( M ) + ..., 



(4.100) 



(4.101) 



where the supers cript refers to the power of the soft scale Q. The first term does not depend on /t and 
equals Co in Eq. (4.98). The //-dependence of Cq'(//) can be determined by solving the renormalization 
group equation 

d ' 



dp 



r (-i) +r (o) +T (i) 



0. 



(4.102) 



One also needs one additional input parameter, such as e. g. a 
This leads to 



(2) 



mm 



2 4 



in order to fix the integration constant. 



(4.103) 



It is then easy to verify that the scattering amplitude T( ^ + + is //-independent up to terms 
of order Q 2 . Further, the effective range function is given at this order by 

2 



k cot 5 



4vr 1 
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r (0) 
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T (o) 



Ti- 1 ) \ Ti- 1 ) 



+ ik , 



(4.104) 



which can be used to predict the "chiral" expansion for the coefficients in the ERE. Here I list only the 
result for the effective range which is sufficient for our purposes. The expressions for V2,3 can be found 
in 161 
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(4.105) 
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Figure 12: Effective potential and scattering amplitude in the Weinberg-like approach. The dashed- 
dotted line refers to the full long-range interaction. Solid dot and filled rectangle refer to the leading 
and subleading contact interactions, respectively. For remaining notation see Fig. 1. 



As expected, the first three terms in the "crural" expansion of r are correctly reproduced at N 2 LO being 
protected by the LETs. The same holds true for the shape parameters Vi, see Ref. |67j . The knowledge 
of atx] for one particular Xj is sufficient to predict ai\ for all k ^ j. 

4.5.3 Weinberg-like approach with a finite cutoff 

An elegant EFT formulation like the one described above which respects the manifest power counting 
at every stage of the calculation is not available in the realistic case of nucleon-nucleon interaction. 
Here, one lacks a regularization prescription for all divergent integrals resulting from iterations of the 
potential in the LS equation which would keep regularization artefacts small without, at the same 
time, introducing a new hard scale in the problem. Contrary to the considered model, the l-7r-exchange 
potential is non-separable and cannot be analytically resummed in the LS equation. In the context of 
chiral EFT for few-nucleon systems, the divergent integrals are usually dealt with by introducing an UV 
cutoff A, which needs to be taken of the order A ~ m s or higher in order to keep regularization artefacts 
small. Clearly, cutoff-regularized diagrams will not obey dimensional power counting anymore. This, 
however, does not mean a breakdown of EFT since power counting is only required for the renormalized 
amplitude. I now consider the Weinberg-like formulation in which the effective potential, given by the 
long-range interaction and a series of contact terms, is iterated in the LS equation to all orders, see the 



work by Lepage [73] for a related discussion. This is visualized in Fig. 12 I carry out renormalization by 
literally following the steps outlined in Ref. [73] and summarized in Ref. [75] in the following way: "The 
theory is fully specified by the values of the bare constants ... once a suitable regularization procedure 
is chosen. In principle, the renormalization program is straightforward: one calculates quantities of 
physical interest in terms of the bare parameters at given, large value of (ultraviolet cutoff) A. Once a 
sufficient number of physical quantities have been determined as functions of the bare parameters one 
inverts the result and expresses the bare parameters in terms of physical quantities, always working at 
some given, large value of A. Finally, one uses these expressions to eliminate the bare parameters in all 
other quantities of physical interest. Renormalizability guarantees that this operation at the same time 
also eliminates the cutoff. " Notice that by iterating the truncated expansion for the effective potential 
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in the LS equation one unavoidably generates higher-order contributions without being able to absorb 
all arising divergences into redefinition of the LECs present in the considered truncated potential. 
Thus, for the case at hand, cutoff dependence in observables is expected to be eliminated only up to 
the considered order in the EFT expansion. I further emphasize that expressing the bare parameters 
(i. e. LECs Cj) in terms of physical quantities is a non-trivial step as the resulting nonlinear equations 
for {Cj} do not necessarily possess real solutions for all values of A, especially when it is chosen to be 
considerably larger than m s . 

To be specific, consider the effective potential at subleading order in the Weinberg- like approach as 



depicted in Fig. 12 



V eS (p\ p) = vi Fi(p') Fi(p) + Co . (4.106) 
The off-shell T-matrix T(p', p; k) can be easily calculated by solving the 2x2 matrix equation 

t(k) = v cS + v cS G(k)t(k) (4.107) 

where I have defined 

Veffip', p) = 7 V) ^eff lip), Tip', p, k) = 7 T (p') tik) 7 (p) , (4.108) 

with 

The integral Ii(k) is given by 
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and is ultraviolet-finite. The divergent integrals Ji(fe) and In(k) are regularized by means of a finite 
cutoff A: 
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2(A; 2 + mf) V2A/ 2 (fc 2 + mf) J y 1 ' V ; 

where s = (^2^Jm 2 — m 2 /m s ^j arccot (mi/ \J m 2 — mf ^ . Neglecting, for the sake of simplicity, the 

finite cutoff artefacts represented by the C(A~ 1 )-terms in Eq. (4.110) and performing straightforward 
calculations, one obtains for the scattering length: 

Trm s {C m [2ai im s (A - smi) + 2m 2 ln(m s /2A)) + vrm/m s ] + 47r 2 a/m s } 

a A = ? ^ • (4.111) 

mi < 2-nm 2 . (ComA + 2-/r 2 ) — Comm/a/ [sm s — 2m/ ln(m s /2A)] > 
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Renormalization is carried out by matching the above expression to the value of the scattering length 
in the underlying model which is regarded as data, 



^underlying 



mi (2a; -l)a s - onm s 
mi (miaia s - m s ) 



(4.112) 



and expressing Co (A) in terms of a un deriying- A straightforward calculation yields the following renor- 
malized expression for the effective range: 
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(4.113) 



In agreement with the LETs discussed above, one observes that the subleading terms in the "crural" 
expansion of r (and Vi, see [67J) are correctly reproduced once Co is appropriately tuned. Notice 
that the smallness of the subleading correction to r\ due to the Co-term in the effective potential as 



compared to the leading contribution given by the first term on the right-hand side of Eq. (4.113) 
is only guaranteed after carrying out renormalization by properly tuning Co (A). The sub-subleading 
and higher-order terms in the "chiral" expansion of r and Vi are not reproduced correctly being not 
protected by the LETs at the considered order. Moreover, since the included LEC is insufficient to 
absorb all divergencies arising from iterations of the LS equation, nothing prevents the appear ance of 



(>2) 



The results in Eq. (4.113) 



positive powers or logarithms of the cutoff A in the expressions for a 
show that this is indeed the case. The dependence on A occurs, however, only in contributions beyond 
the accuracy of calculation and, obviously, does not affect the predictive power of the EFT as long as 
the cutoff is chosen to be of the order of the characteristic hard scale in the problem, A ~ m s . 

An important misconception that appears frequently in the literature is related to the treatment of 
the cutoff by employing very large values of A or even regarding A — > oo. While this is perfectly fine 
in ChPT, where observables are calculated perturbatively and all emerging UV divergencies can be 
absorbed by the corresponding counterterms at any fixed order in the chiral expansion, this is not a 
valid procedure for the case at hand. Let us further elaborate on this issue using the above example. 
At first sight, the appearance of positive powers of A and/or logarithmic terms in the predicted "chiral" 



expansion of the subthreshold parameters, see Eq. (4.113), may give a (wrong) impression that no finite 
limit exists for r\ and (vi)\ as A — > oo. Actually, taking the limit A — > oo does not commute with the 
Taylor expansion of the ERE coefficients in powers of m;. Substituting the value for Co (A) resulting 



from matching Eq. (4.111) to (4.112) into the solution of the LS equation (4.107) and taking the limit 
A — > oo yields the following finite, cutoff-independent result for the inverse amplitude: 
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(4.114) 
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The corresponding infinite-cutoff prediction for the effective range has the form: 
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(4.115) 

where the ellipses refer to O (mf^j -terms. One observes that the result after removing the cutoff fails to 

reproduce the low-energy theorem by yielding a wrong value for a?' . This also holds true for the a^) 
|67j . Notice that, by construction, the scattering length is still correctly reproduced. The breakdown of 
LETs in the Weinberg-like approach in the A — > oo limit can be traced back to spurious A-dependent 
contributions still appearing in renormalized expressions for observables, see Eq. (4.113), which are 
irrelevant at the order the calculations are performed in the regime A ~ m s but become numerically 
dominant if A ^> m s . Due to non-renormalizability of the effective potential as discussed above, such 
spurious terms do, in general, involve logarithms and positive powers of A which, as A gets increased 
beyond the hard scale m s , become, at some point, comparable in size with lower-order terms in the 
"chiral" expansion. For example, the appearance of terms linear in A would suggest the breakdown of 
LETs as the cutoff approaches the scale A ~ m? s / mi . The unavoidable appearance of ever higher power- 
law divergences when going to higher orders in the EFT expansion implies that the cutoff should not be 
increased beyond the pertinent hard scale in Weinberg-like or Lepage-like approach to NN scattering 
leading to A ~ m s as the optimal choice. It is furthermore instructive to compare the predictions for 
the effective range in Eqs. (4.105) and (4.113) corresponding to two different renormalization schemes. 
One observes that taking A 3> m s in Eq. (4.113) has an effect which is qualitatively similar to choosing 
fi> mi in Eq. (4.105) and corresponds to an improper choice of renormalization conditions in the EFT 
framework. 



4.5.4 Toy model with local interactions: a numerical example 

Having clarified the important conceptual issues related to nonperturbative renormalization in the 
context of chiral EFT for two nucleons, I now turn to the last toy-model example and give some 
numerical results. 

Consider two nucleons interacting via the local force given by a superposition of two Yukawa potentials 
corresponding to the (static) exchange of the scalar light and heavy mesons of masses mi and m s , 
respectively: 

-m ; r p-m a r 

V( r ) = ai—— + a s (4.116) 

This type of potentials is sometimes referred to as Malfliet-Tjon potential. Motivated by the realistic 
case of the two-nucleon force, I choose the meson masses to be mi = 200 MeV and m s = 750 MeV. 
Further, I adjust the dimensionless strengths ati s in such a way that the potential features an S-wave 
bound state ("deuteron") with the binding energy Eb = 2.2229 MeV. A suitable combination is given 
by oli = —1.50 and a s = 10.81. With the parameters specified in this way, the potential is depicted in 
Fig. [T3| The corresponding momentum-space potential can be easily obtained by making the Fourier 
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Figure 13: Toy-model potential in Eq. (4.116). The dashed (dashed-dotted) line depicts the short-range 
(long-range) part proportional to a s (aj) while the full potential is shown by the solid line. 



transformation: 

V(p , p ) = 2 + ^2~, 2 • (4.117 

q 2 +mf q z +mj 

Here, q = p' — p denotes the momentum transfer. I treat the nucleons in this example as identical, spin- 
less particles. Thus, in the partial wave basis, the only nonvanishing matrix elements {l',j,p' \V\l,j,p) 
correspond to I = I' = j. I only consider the S-wave here. Because of no spin dependence, the matrix 
element Vq(j)', p) = (0,0, p' | V|0, 0, p) can be obtained by simply integrating over the angle 8 between 
p' and p : 

r+l 

V (p',p) = 2nJ d(case)V(p',p,6) 

4vr 2 . f(p' +p) 2 + mf\ 4?r 2 n p ' +p f + m ^ 



= on— In — rs- 2 )+a s — In, 

p p V \P — p) + m t J p p V \P ~ p) + m s 

= VZ(p',p) + V s (p',p). (4.118) 
Contrary to the previously considered case of a separable interaction, the LS equation 

r l 2 dl m 
T (p, p'; k) = V (p, p') + J (^3^o(p, l) k2 _ l2 + ie T o(l, P'; k) , (4.119) 

cannot be solved analytically for the Malfhet-Tjon-type potentials. It can, however, be solved numeri- 
cally using the standard methods, see e.g. |12j . With the parameters specified above, one obtains the 
phase shift which is shown by the solid line in the left panel of Fig. 14 It is fairly similar to the 
neutron-proton 3 Si phase shift, cf. the left panel of Fig. [lj 
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Figure 14: Left panel: phase shifts resulting from the original and effective potentials in Eqs. (4.117) 



and (4.120), respectively. Right panel: Lepage plot showing the absolute error in phase shifts at various 



orders in the low- momentum expansion versus lab. energy. 



I now develop an effective potential that describes the same physics as the underlying one at momenta 
of the order of Q ~ m/. Up to N 2 LO, it takes the formj^] 



Vf{ P \ p) = V l (p', p) + [Co + C 2 ( P ' 2 + P 2 ) + C 4P ' 2 p 2 ] f A (p', P ) 



(4.120) 



The regulator function prevents the appearance of ultraviolet divergences in the LS equation and is 
chosen of the form: 



fA(p',p) = exp 



p'2 _|_ p2 
A 2 " 



(4.121) 



I set the cutoff A = 500 MeV, solve the LS equation with the effective potential V cS (p', p) and adjust 
the LECs to reproduce the coefficients in the ERE as follows: 

• LO: C2 = C4 = 0, Co is tuned to reproduce a; 

• NLO: C4 = 0, Co,2 are tuned to reproduce {a, r}; 

• N 2 LO: Co,2,4 are tuned to reproduce {a, r, V2}- 



With the LECs being fixed as described above, the predictions for the S-wave phase shift at various 
orders in the low- momentum expansion are summarized in the left panel of Fig. |14[ The so-called Lepage 
plot in the right panel of this figure shows absolute errors in the phase shift, A5(E\ ai y ) ) = <J U nder lying - ^effj 
versus energy. It is plotted in radians. One reads off from this plot that the laboratory energy, at which 
the expansion breaks down, is of the order of £i a b = 2k 2 /m ~ 250 MeV. This corresponds to the 

11 One can write down another contact interaction with four derivatives whose matrix elements in the on-shell kinematics, 
i.e. with p — p' , cannot be disentangled from the C4-term. 
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momentum scale of the order of A ~ 350 MeV, in a good agreement with the expected breakdown 



scale of the modified effective range expansion of the order of m s /2, see section 4.2 Notice that the 



effective theory is, as desired, able to go beyond the ERE, whose range of convergence is indicated by 



the vertical lines in Fig. 14 The "deuteron" binding energy is found to be reproduced correctly with 
5 significant digits at N 2 LO: 

E%° + 6E% LO + 5E^ LO = 2.1594 + 0.638 - 0.0003 = 2.2229 MeV . (4.122) 

For further illustrative quantum mechanical examples and a discussion on renormalization in the context 
of the Schrodinger equation, the reader is referred to the excellent lecture notes by Lepage [74"] . 

5 Nuclear forces from chiral EFT 

In this section I outline and exemplify some methods which can be used to derive nuclear forces from 
chiral EFT. 

5.1 Derivation of nuclear potentials from field theory 

The derivation of a potential from field theory is an extensively studied problem in nuclear physics. 
Historically, the important conceptual achievements in this field have been done in the fifties of the 
last century in the context of the so-called meson field theory. The problem can be formulated in the 
following way: given a field theoretical Lagrangian for interacting mesons and nucleons, how can one 
reduce the (infinite dimensional) equation of motion for mesons and nucleons to an effective Schrodinger 
equation for nucleonic degrees of freedom, which can be solved by standard methods? It goes beyond 
the scope of this work to address the whole variety of different techniques which have been developed 
to construct effective interactions, see Ref. [76] for a comprehensive review. I will now briefly outline a 
few methods which have been used in the context of chiral EFT. Similar methods are frequently used 
in computational nuclear physics in order to reduce a problem to a smaller model space which can be 
treated numerically. 

I begin with the approach developed by Tamm [TTJ and Dancoff |78| which in the following will be 
referred to as the Tamm-Dancoff (TD) method. Consider the time-independent Schrodinger equation 

(H + Hi)\*)=E\V), (5.123) 

where \^>) denotes an eigenstate of the Hamiltonian H with the eigenvalue E. One can divide the full 
Fock space in to the nucleonic subspace \ <p) and the complementary one \tp) and rewrite the Schrodinger 



equation (|5.123|) as 

r/Hr] r/HX 
XHi] XHX I V |$ 



(5.124) 
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where I introduced the projection operators r\ and A such that \<p) = rf\^f), \ip) = A|\l/). Expressing the 
state \ip) from the second line of the matrix equation (5.124) as 

H\4>) , (5.125) 



E-XHX 



and substituting this in to the first line, one obtains the Schroedinger-like equation for the projected 
state \4>): 

(# o + Kff D (£))l0> = £|0>, (5-126) 
with an effective potential V e s(E) given by 

V£ D (E) = r,H ir , + vHi\ E _ 1 XHX XH lV . (5.127) 



This definition of the effective potential corresponds exactly to the one given in section 4.1 in the 
context of "old-fashioned" time-ordered perturbation theory. To evaluate V^P(E) one usually relies 
on perturbation theory. For example, for the Yukawa theory with a single ttNN vertex Hj = gH\, 
VtfP(E) up to the fourth order in the coupling constant g is given by 



K£ D (£) = -v' 



g 2 H 1 ^^H l +g 4 H 1 ^^H l ^^H 1 ^^H 1 + 0{g Q 

» 1 tt TP ^ TJ TP TJ TP TJ TP Ki) 

Hl() — Hi HtQ — Hi I1Q — Hi Hlfl — Hi 



rj, (5.128) 



where the superscripts of A refer to the number of mesons in the corresponding state. It is important to 
realize that the effective potential V c q(E) in this scheme depends explicitly on the energy, which makes 
it inconvenient for practical applications (especially for calculations beyond the two-nucleon system). 
In addition, the projected nucleon states \(f>) have a different normalization compared to the states \^S?) 
we have started from (which are assumed to span a complete and orthonormal set in the whole Fock 
space) 

(&|^-) = (tfiltfj) - WMj) = Sij - (falHjX ( g _ 1 AffA ) , (5.129) 

since the components do, in general, not vanish. 

The above mentioned deficiencies are naturally avoided in the method of unitary transformation |79} [80] . 
In this approach, the decoupling of the rj- and A-subspaces of the Fock space is achieved via a unitary 
transformation U 

H^U^HU=( f*" .1 V (5.130) 



XHX 

Following Okubo |79j . the unitary operator U can be parametrized as 

/ 1/(1 + A* A)- 1 /* -A* (I + AAt)~W \ 
u - { 4(i + AU)-i/J A(l + ^t)-V2 J , 

with the operator A = XArj. The operator A has to satisfy the decoupling equation 

X(H — [A, H] — AHA) 77 = (5.132) 
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in order for the transformed Hamiltonian H to be of block-diagonal form. The effective 77-space potential 
V^g T can be expressed in terms of the operator A as: 



KfT T = V{H 



(1 + rfA)-V 2 {H + A*H + HA + A ] HA){1 + A ] A)' 1 / 2 - H 



(5.133) 



The solution of the decoupling equation and the calculation of the effective potential according to 



Eq. (5.133) can be carried out perturbatively in the weak-coupling case. For the previously considered 



case of the Yukawa theory, the decoupling equation can be solved recursively by making the ansatz 

oo 

A = ^g n A i - n \ (5.134) 



n=l 



The resulting effective potential takes the form: 



2 / 

-9 V 



1 \i 

-Hi— —Hi + h. c. 

2 Hq — E„ 



4 / 

rj-g 7? 



l Hi /TT X \ . HtfjHv 



A 1 



2 (Hq — E n i 



A 1 



8 (Hq — E v > 



-Hi 



-HifjHi- 



Ho — Efj)(Ho — E v 
A 1 



-Hi^—Hi 

2 (Hq — E v ) (Hq — E v ) (Hq — E v 
Hi 



Ho — Efj)(Ho — E v 
A 1 



(Hq — E v ')(Hq — Ef f 



-HifjHi- 



Hq — Efj 



Hi 

-Hi + h. c. 



V + 0(g 6 ). 



(5.135) 



Here, rj, rj and fj denote projection operators onto the purely nucleonic states. Different notation is 
only used to indicate what state the energies in the denominators correspond to. In contrast to V^ D , 
V^j T does not depend on the energy E which enters the Schrodinger equation. Another difference to the 
Tamm-Dancoff method is given by the presence of terms with the projection operator fj which give rise 
to purely nucleonic intermediate states. These terms are responsible for the proper normalization of 
the few-nucleon states. In spite of the presence of the purely nucleonic intermediate states, such terms 
are not generated through the iteration of the dynamical equation and are truly irreducible. Since all 
energy denominators entering V^g T correspond to intermediate states with at least one pion, there is 
no enhancement by large factors of m/Q that occurs for reducible contributions. 

Exercise: 

1. Calculate A 1 , A^ and A^ by solving the decoupling equation for the considered case of Yukawa 



theory and verify the expression for V^ T in Eq. (5.135). 



2. Consider the disconnected Feynman diagram in Fig. 



15 and draw all possible time-ordered diagrams. 



Using Eqs. (5.128) and (5.135) show that, in contrast to the TD approach, these diagrams do not 
contribute to the nucleon-nucleon potential in the method of unitary transformation. Use the static 



approximation for the nucleons in order to simplify the calculations (i.e. set: E = E v = E„ 



Efj 



0). 



47 



Figure 15: An example of a disconnected diagram that does not contribute to the NN potential in the 
method of unitary transformation. 



The two methods of deriving effective nuclear potentials are quite general and can, in principle, be 
applied to any field theoretical meson-nucleon Lagrangian. In the weak-coupling case, the potential can 
be obtained straightforwardly via the expansion in powers of the corresponding coupling constant (s). 
For practical applications, it is helpful to use time-ordered diagrams to visualize the contributions to 
the potential, see Fig. |6j In "old-fashioned" perturbation theory or, equivalently, the Tamm-Dancoff 
approach, only irreducible diagrams contribute to the potential. In the method of unitary transforma- 
tion one can draw both irreducible and reducible graphs whose meaning, however, differs from that of 
diagrams emerging in time-ordered perturbation theory. The coefficients in front of various operators 
and the energy denominators can, in general, not be guessed by looking at a given diagram and have 



to be determined by solving the decoupling equation (5.132) for the operator A and using Eq. (5.133). 



Application of the above methods to the effective chiral Lagrangian requires the expansion in powers 
of the coupling constants to be replaced by the chiral expansion in powers of Q/A x . This issue will be 
dealt with in the next section. 



5.2 Method of unitary transformation 



To apply the method of unitary transformation to derive nuclear forces in chiral EFT it is useful to 



rewrite the power counting discussed in section 4.1 into a different form which is more suitable to carry 
out algebraic manipulations described above. 

We begin with Weinberg's original power counting expression for A-nucleon diagrams involving C 
separately connected pieces: 



4 _ N + 2 (L - C) + V ^ > A i 



di + -rii 



2. 



(5.136) 



This expression is a generalization of Eq. (4.58) to the case C > 1. Its derivation can be found in 
Ref. 



[81j. There is one subtlety here that needs to be addressed: according to Eq. (5.136), the chiral 
dimension v for a given process depends on the total number of nucleons in the system. For example, 
one-pion exchange in the two-nucleon system corresponds to N = 2, L = 0, C = 1 and T^Aj = 
and, therefore, contributes at order v = 0. On the other hand, the same process in the presence of a 



48 



third (spectator) nucleon leads, according to Eq. (5.136), to v = —3 since N = 3 and C = 2. The origin 
of this seeming discrepancy is due to the different normalization of the 2N and 3N states: 



S 3 (pi -pi )5 3 {p 2 ' ~P2 



2N : (pi^IpYpY) 
3N : {piP2Pz\pi P2 Pz) = $ 3 (pi -pi )8 3 (p2 -P2)5 3 (p3 -p 3 ) 



(5.137) 



It can be circumvented by assigning a chiral dimension to the transition operator rather than to its 
matrix elements in the N- nucleon space. Adding the factor 3iV — 6 to the right-hand side of Eq. ( 5.136 ) 



in order to account for the normalization of the N- nucleon states and to ensure that the LO contribution 
to the nuclear force appears at order v = we obtain 



-2 + 2N + 2(L - C) + v i A i 



(5.138) 



This expression provides a natural qualitative explanation of the observed hierarchy of nuclear forces 
V2N > V3N 3> V4N • • • with 



V2N 
V3N 
V4N 



y(0) , y(2) , F (3) , ta(4) 
l/ 2N + V 2N + K 2N + l/ 2N + 

> / 4N + • • • 1 



(5.139) 



as shown in Fig. 16 



The form of power counting in Eq. (5.138) is still of less use for our purpose since the resulting chiral 



dimension is given it terms of the topological quantities such as N, C and L which is not appropriate for 
algebraic approaches such as the method of unitary transformation. Using certain topological identities, 



see [82], Eq. (5.138) can be rewritten in a more suitable form: 



-2 + J>* 



h . 



di + -Hi + pi 



4. 



(5.140) 



The quantity Ki which enters this expression is nothing but the canonical field dimension of a vertex of 
type i (up to the additional constant —4) and gives the inverse mass dimension of the corresponding 
coupling constant. In fact, this result can be obtained immediately by counting inverse powers of the 
hard scale A x rather than powers of the soft scale Q (which is, of course, completely equivalent). Indeed, 
since the only way for the hard scale to be generated is through the physics behind the LECs, the power 
v is just the negative of the overall mass dimension of all LECs. The additional factor —2 in Eq. (5.140) 
is a convention to ensure that the contributions to the nuclear force start at v = 0. I encourage the 



reader to verify the equivalence of Eqs. (5.140) and (5.138) for specific diagrams. One immediately 



reads off from Eq. (5.140) that in order for perturbation theory to work, the effective Lagrangian must 
contain no renormalizable and super-renormalizable interactions with Kj = and Kj < 0, respectively, 
since otherwise adding new vertices would not increase or even lower the chiral dimension v. This 
feature is guaranteed by the spontaneously broken chiral symmetry of QCD which ensures that only 
non-renormalizable interactions enter the effective Lagrangian. 
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Figure 16: Diagrams that give rise to nuclear forces in ChEFT based on Weinberg's power counting. 
Solid and dashed lines denote nucleons and pions, respectively. Solid dotes, filled circles and filled 
squares and crossed squares refer to vertices with Aj = 0, 1, 2 and 4, respectively. 



While Eq. (5.140) does not say much about the topology and is, therefore, not particularly useful to 
deal with diagrams, it is very convenient for algebraical calculations. In fact, it formally reduces the 
chiral expansion to the expansion in powers of the coupling constant, whose role is now played by the 
ratio Q/A. Applying the canonical transformation to the chiral Lagrangian and writing the resulting 
Hamiltonian in the form 



(5.141) 



the operator A can be calculated by solving Eq. (5.130) recursively, 

oo 



A 



A (a) 



(5.142) 



a=l 



1 



Erj — E\ 



A 



a-1 



a-1 



a-2 a—j—1 

i=l i=l i=l j=l 



The expressions for the unitary operator and the effective potential then follow immediately by substi 
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Figure 17: Diagrams that correspond to the operator in Eq. (5.144). Graph (a) yields a single-body 
contribution to the nuclear Hamilton operator while diagrams (b) and (c) give rise to the l-7r-exchange 
NN potential. 



tuting Eqs. (5.141) and (5.142) into Eq. (5.133). 



5.3 The lit- and the leading 27r-exchange potentials 



I now illustrate how the above ideas can be applied in practice. I begin with the simple case of the 
l7r-exchange potential at leading order, i.e. v = 0. The only relevant contribution to the interaction 
Hamilton density is given by 



9A 
2F n 



N^a ■ (V7r • t)N . 



(5.143) 



where the superscript of T~L gives the canonical dimension Ki defined in Eq. (5.140). The relevant 



operator that contributes to the effective Hamiltonian after performing the unitary transformation is 



given by the first two terms in Eq. (5.135): 



yVT 



(5.144) 



where ui denotes the pion free energy and I made use of the static approximation as appropriate at 
LOj^] Notice that V^ T in the above equation also gives rise to a one-body operator that contributes 
to the nucleon mass shift, see graph (a) in Fig. 17 To compute the expression for the l-7r-exchange 



potential I first express the pion and nucleon fields in the interaction picture in terms of the creation 
and destruction operators: 



N(x) 



e- l «- x ai (k) + e lk - x a 



f d 6 k 1 r 
J (2vr)3/2 ^ET k [ e 



(5.145) 



Corrections to the static l7r-exchange potential are suppressed by Q 2 /m 2 
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where ojk = yk 2 + m^ and v (e) denotes a Pauli spinor (isospinor). ai(k) and a\(k) denote a de- 
struction and creation operator of a pion with isospin i. Further, bt(p, s) (b\(p, s)) is the destruction 
(creation) operator of a non-relativistic nucleon (i.e. po = p 2 /(2m)) with the spin and isospin quantum 
numbers s and t and momentum p. The creation and destruction operators of the pion (nucleon) field 
satisfy the usual commutation (anti-commutation) relations. The l7r-exchange potential can be calcu- 



lated by substituting the expressions in Eq. (5.145) for pion and nucleon fields into Eq. (5.144) and 
evaluating the matrix element 



{a' x a' 2 \ 



OJ 



V |aicn2) 



1 



(5.146) 



Here, P (P') denotes the total momentum of the nucleons before (after) the interaction takes place. 
Further, a% and ol i denote collectively the initial and final quantum numbers of the nucleon i (momen- 
tum, spin and isospin). To keep the expressions for the potential compact, they are commonly given in 
the operator form with respect to the spin and isospin quantum numbers using the corresponding Pauli 
matrices at and Tj of a nucleon i. A straightforward calculation yields the final form of the l-7r-exchange 
potential: 



V 2N 



9a <h ■ gjj 
4P2 



.. ,- + Mr ri - T2 - (5 - 147) 

Clearly, this familiar result for the static l-zr-exchange potential can be obtained in a much simpler 
way by evaluating the corresponding Feynman diagram since it does not generate reducible topologies. 
One-loop corrections to the static l-7r-exchange potential and renormalization within the method of 
unitary transformation are discussed in detail in Ref. |83j . Notice further that when calculating the 



matrix element in Eq. (5.146), I discarded the contributions corresponding to graph (a) in Fig. 5.144 
with one of the nucleons being a spectator and the contributions from diagrams (b) and (c) with the 
nucleon labels a' x and a' 2 being interchanged. The latter emerge automatically when the potential is 
inserted into the corresponding dynamical equation due to the antisymmetric nature of the two-nucleon 
wave function. 

As a final example, I discuss the leading 2-7r-exchange potential arising at order v = 2 from the box and 



crossed box diagrams (the last two diagrams in the second raw in Fig. 16 Again, the only vertex we 



need is given in Eq. (5.143). The relevant operators that contribute to the effective nuclear Hamiltonian 



after performing the unitary transformation are listed in Eq. (5.135) 



UT 



A 2 



cff 



OJ 0J\ + 0J2 



H 



(i) 



oj 2 



oj- 



OJ 



+ 



1^(1)^(1)^(1)^(1) 
2 oj oj 1 



(5.148) 



The contribution to the 27r-exchange potential results from evaluating the matrix element (a^a^l ^ ! ff T | a i Q; 2) 
which can be computed along the same lines as above. Calculations of that kind can be optimized by 
using a diagrammatic approach and formulating a sort of "Feynman" rules. The building blocks are 
given by vertices and energy denominators that play the role of propagators in Feynman diagrams. 
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Consider, for example, time-ordered box diagrams (b)-(g) in Fig. [6j All these graphs have an identical 
sequence of non-commuting vertices generating exactly the same isospin-spin-momentum structure in 
the resulting potential. Thus, the energy denominators for different diagrams arising from the operators 



in Eq. (5.148) can be added together yielding the result 



} uj + UJiUJ 2 + u>\ 



(5.149) 



The same result but with an opposite sign is obtained for the sum of the energy denominators for the 
crossed-box diagrams. 



Exercise: show that the operators in Eq. (5.148) do not give rise to the 2ir-exchange three-nucleon force. 



What would be the result for the three-nucleon force if one would employ time- ordered perturbation 
theory (in the static approximation) instead of the method of unitary transformation? 



The vertex in Eq. (5.143) gives rise to the "Feynman" rule 

1 



f 9 A 
''2F n 



/ 2oj n 



(5.150) 



for an incoming (outgoing) pion with momentum q (— q) and the isospin quantum number a. Here, i is 
the nucleon label. Putting everything together, we end up with the contribution from the box diagram 
of the form 



T r2n, box / ->\ 



(2tt) j 



9A 
2F n 



T\ ■ T2T1 ■ T2 



d 3 h d% _ r _ r . f _ f 
(2vr)3(2vr)3 °1 • W W • 



1 1 OJ? + OJU 0Jl o + uf 



(5.151) 



where the factor (2ir) 3 in front of the integral is due to the normalization of the potential, see Eq. (5.146). 
The contribution of the crossed-box diagrams can be written as 



V, 



2-7T, cr.— box / -»\ 
2N (?) 



-(27T) 



_9A_ 

2 Fir 



E„ / d 3 h d 3 l 2 _ - 



(2vr) 3 (2vr) 3 



1 1 2 < + ^ + < 8(l 1 + l2 



2u h 2u> h u^uf^u^ + ui 2 ) 



(5.152) 



Adding the two expressions together and performing straightforward simplifications one obtains the 
total contribution to the leading 2-7r-exchange proportional to g\: 



t r2-7r, total / ->\ 



dH 



32F*J (2tt) 



A 



n • r 2 [I 



+ 6ai-qxla2-qxl 



(5.153) 
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with uj± = y (q±l) 2 + 4M^. The integrals appearing in the above expressions are ultraviolet divergent 
and need to be regularized. This can be achieved using standard methods such as e.g. dimensional 
regularization. Cutoff regularization can be applied equally well. In the infinite-cutoff limit, A — > oo, 
the regularized integrals can be decomposed into a finite non-polynomial part (with respect to the 
momentum transfer q ) and polynomial in momenta terms that may diverge as A goes to infinity. Such 
a decomposition follows from the local nature of the ultraviolet divergences and implies the uniqueness 
of the non-polynomial part (in the limit A — > oo). This makes perfect sense from the physics point of 
view since the nonpolynomial part of the potential controls its long-range behavior which should not 
depend on the details of regularization at short distances. For the non-polynomial parts of the relevant 
integrals one obtains: 

I fZ* , ^,\„ S = ^(4M^ + , 2 )L(,) + ..., 

1 512M^ + 384M*q 2 + 156M 2 q 4 + 23q e 



dH 


f 2 


(2tt) 3 


oj + ijJ-(uj + + o;_) 


dH 


f 4 + q 4 


(2tt) 3 




d 3 l 




(2vr) 3 





60vr 2 4M2 + q 2 L ^) + • • • ' 

J 3= I 7^ n = +••-. ( 5 - 154 ) 



where q = \q\ and the ellipses refer to terms polynomial in q. Note that ^3 does not give rise to any 
non-polynomial terms. Further, I have introduced the loop function L(q) defined as: 

tM^l. ^f . (5.155) 



Using the identity 



uii. + uj+uj- + u)_ Id 



u\urt{u+ +u-) 2 d(M 2 ) uj + uj-{u] + +uj-) 



(5.156) 



one can express all integrals entering Eq. (5.153) in terms of Ji,2,3 as follows 
d 3 l uPl + oj + uj^ + uPl /t9 2 .2 1 d , 2 



( l2 -1 2 ) = -o^^-V/i), (5-157) 



(2tt) 3 w 3 w 3 (w+ + w_) v h > 2 d(M 2 ) 

dH U% + U + LU_ +^ 2 ;; 1 d r / 1 r \ . / 1 3, 



(2tt) 3 w 3 wi(w + +a;_) tJ 4 d(M|) U q* a ) %3 \q 

The final result for the 27r-exchange potential oc g\ then takes the form: 

- - 3^ - ' - ( 20M - 2 + 23 « 2 + ilr?) £ « 

3o 4 

- • ^2 • q~ q 2 #1 • B 2 ) L(q) + .... (5.158) 
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The polynomial in momenta, divergent (in the limit A — > oo) terms have the form of contact interactions 
that are anyway present in the potential at a given order and can be simply absorbed into an appropriate 
redefinition of the LECs C,. 

Exercise: verify the result for the non-polynomial part of V 2 2 ^' tot&l using dimensional regularization. 
Use the equality 

1 2 f 00 1 1 

7 \=-/ d P~2 7y?—2 59. ( 5 - 159 ) 

CJ + U)-(u> + + vr J r UJ 2 _ + f3 2 u\ + f3 2 

to get rid of the square roots in the integrand. The resulting integrals can be dealt with in the usual 
way by introducing the corresponding Feynman parameters. 

In coordinate space, contact interactions have the form of the delta function at the origin, 5(f), and 
derivatives thereof. In contrast, the nonpolynomial pieces give rise to the potential at finite distances. 
To see this let us take a closer look at the obtained expression for the 2-7r-exchange potential. First, it 
should be emphasized that the Fourier transformation of the nonpolynomial terms alone is ill defined 
since they grow as q goes to infinity. The potential V^ ,tot (f) at a finite distance, r / 0, can be 
obtained from V 2 2 ^ ,total (q ) via 

V 2 Y° ta \r) = J (0 e~^V 2 Y 0ta \q) F (|) , (5.160) 

where F(q/A) is an appropriately chosen regulator function such as e.g. F = exp(— q 2 /A 2 ). Alterna- 
tively, one can use a (twice-subtracted) dispersive representation by expressing the potential V 2 ^ totai (q) 
in terms of a continuous superposition of Yukawa functions. For example, the central part of the po- 
tential in Eq. (5.158) can be written as (HU [85] 



t .2tt, central, s f°° , 1 P(lA - ...... 

V 2N W = / d ^ 2 _l 2 ' (5.161 



where the spectral function p{q) is given by 

p(M) = Im[^ ccntral (0+-^) 

9 a ( 2 2 48M^ \ V V - AM 2 

The Fourier transformation can be easily carried out in this spectral representation by first integrating 
over q and then over the spectrum fi. This leads to the central potential of the form 

V 2 - N ^\r) = ~ 12 g^4 r 4 ti • T2 [(23 + 12x 2 )K 1 (2x) + x(23 + ix 2 )K (2x)] , (5.163) 

where Ki denote the modified Bessel functions and x = M n r. At large distances, the potential behaves 
as exp( -2M n r)/r 3 / 2 . The expressions for the remaining components of the 2-7r-exchange potential up to 
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Figure 18: A schematic picture of the two-nucleon potential derived in chiral EFT in a given partial 
wave. 



the chiral order Q 3 , both in momentum and coordinate space, can be found in Ref. |84j . The order-Q 4 
contributions are given in Ref. |86j . The expressions for pion exchange potentials derived in chiral 
EFT at large distances are controlled by low values of n for which the chiral expansion is expected 
to converge. At shorter distances, the large-/x components in the spectrum start to contribute which 
cannot be computed reliably in chiral EFT. This is visualized in Fig. 18 An extended discussion on 
the resulting theoretical uncertainty can be found in Ref. [85J. 



It is instructive to compare the toy models considered in section 4.5 with the nucleon-nucleon potential 



derived in chiral EFT whose structure is symbolically illustrated in Fig. 18 The main conceptual 



difference is due to the lack of an exact (regular) expression for the long-range force in the realistic 
case of nucleon-nucleon interaction. Rather, it is represented in terms of the chiral expansion of the 
pion-exchange potential which is valid at large distances and behaves singular as r — > 0. This raised 
debates on the relative importance of the long- and short-range components in the potential and the 
most efficient way to organize the expansion for low-energy observables. There is little consensus on 
this issue in the literature (yet). 

Chiral 27r-exchange potential is, clearly, the most interesting new ingredient of the two-nucleon force 
from the chiral EFT point of view: it is the next-longest-range contribution after the well established 
l7r-exchange potential whose form is strongly constrained due to the chiral symmetry of QCD. Notice 
that three-pion exchange is already considerably less important for low-energy nuclear dynamics. The 
evidence of the chiral 2-7r-exchange potential up to N 2 LO has been verified in the Nijmegen PWA [87J, 
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see Ref. [88J for a similar investigation. In their analysis, the Nijmegen group utilized the long-range 
interaction above some distance b as input in order to constrain the behavior of high partial waves. The 
missing intermediate and short-range components are simulated by suitably chosen energy-dependent 
boundary conditions. The number of parameters entering the boundary conditions needed to achieve 
a perfect description of the data thus may be viewed as a measure of physics that is missing in the 
assumed long-range force. As demonstrated in Ref. |87j . adding the two-pion exchange potential derived 
at N 2 LO in chiral EFT to the l-7r-exchange potential and the appropriate electromagnetic interactions 
allowed for a considerable reduction of parameters (from 31 to 23 for b = 1.4 fm in the case of proton- 
proton scattering) with even a slightly better resulting x 2 . This is a big success of chiral EFT in the 
two-nucleon sector. 

6 Summary 

In these lectures, I have outlined the foundations of chiral effective field theory and the application of 
this theoretical framework to the nuclear force problem. The method allows for a systematic derivation 
of nuclear forces with a direct connection to QCD via its symmetries. These lecture notes are mainly 
focused on the conceptual aspects and do not cover applications of the novel chiral potentials to the 
few-nucleon problem and various related topics such as e.g. isospin breaking effects, few-baryon systems 
with strangeness, electroweak and pionic probes in the nuclear environment, nuclear parity violation 
and chiral extrapolations of few-baryon observables. For a discussion on these and other topics as well 
as for a detailed description of the structure of the two-, three- and four-nucleon forces in chiral EFT 
the reader is referred to recent review articles [89} 190] , see also [M] • There are many frontiers where 
future work is required. These include a better understanding of the power counting in the few-nucleon 
sector, the consistent inclusion of electroweak currents, and the development of chiral EFT with explicit 
A(1232) degrees of freedom. 
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